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ABSTRACT 


THE  EFFECT  OF  SYSTEM  DEPOLARIZATION 
ON  MEASUREMENT  QUALITY 


Robert  Jay  Nelson,  M.S. 

The  University  of  Texas  at  Arlington,  1986 

Supervising  Professor:  Andrew  J.  Blanchard 

Radar  systems  are  susceptible  to  depolarization  effects  caused  by  the  system  and 
the  environment.  System  effects  include  depolarization  caused  by  antennas  and 
measurement  configurations.  Antenna  effects  arise  from  the  assumption  that  the  cross 
polarization  isolation  response  is  infinite  across  the  beamwidth  of  the  antenna.  In  fact,  the 
isolation  response  can  be  as  high  as  -20  dB  at  boresight  to  0  dB  2-3  degrees  off-axis. 
Depolarization  due  to  measurement  configurations  comes  from  the  translation  necessary  to 
match  antenna  coordinate  system  to  the  clutter  coordinate  frame.  For  large  beamwidths  at 
small  angles  of  incidence  this  translation  can  cause  significant  depolarization.  These 
phenomena  have  been  investigated  separately  in  the  monostatic  case. 

This  research  extends  previous  work  done  in  the  monostatic  case  to  the  bistatic 
case.  A  computer  simulation  was  developed  which  models  a  bistatic  clutter  environment. 
The  simulation  models  depolarization  sources  caused  by  system  and  surface  effects.  A 
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linear'y  polarized  wave  is  sent  through  a  transmitting  antenna.  The  coordinate  frame  of  the 


transmitted  wave  is  matched  to  the  clutter  frame.  The  clutter  scatters  a  portion  of  the  wave 
toward  the  receiver.  The  coordinate  frames  are  translated  and  the  wave  proceeds  through 
the  receiving  antenna.  The  result  is  a  "depolarization  figure  of  merit"  attributable  to  the 
system  effects  investigated.  Baseline  runs  were  done  to  establish  the  effect  of  the  individual 
depolarization  source  on  the  degradation  of  system  performance.  Significant  depolarization 
was  found  to  occur  off-axis  due  to  the  individual  contributions  of  both  translation  and 
antenna  effects.  When  combined,  these  effects  can  create  a  large  signal  contribution  located 
off-axis  due  solely  to  system  induced  depolarization.  This  additional  noise  will  serve  to 
hide  any  target  located  in  the  boresight  of  the  antenna.  The  conclusion  reached  from  this 
analysis  is  that  the  inclusion  of  system  induced  depolarization  effects  in  the  modeling  of  a 
radar  performance  is  paramount  in  evaluating  the  system  capability. 
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CHAPTER  ONE 


INTRODUCTION 

The  measurement  of  cross-polarized  scattering  information  is  becoming  more  and 
more  important  in  today's  modem  radar  cross  section  measurement  equipment.  Knowledge 
of  the  scattering  processes,  measurement  geometry,  and  system-induced  depolarization 
effects  is  important  in  order  to  properly  assess  the  quality  of  the  measurements  actually 
made. 

This  thesis  explores  two  different  system  induced  depolarization  effects  to  assess 
their  influence  on  the  quality  of  scattering  measurements  made  in  a  bistatic  environment  A 
bistatic  environment  exists  when  the  transmitting  and  receiving  antennas  are  in  two  different 
locations.  When  the  transmitting  and  receiving  antennas  are  in  the  same  location,  a 
monostatic  environment  exists.  The  two  system  induced  depolarization  effects  studied  are 
translation  depolarization  and  antenna  depolarization.  Translation  depolarization  occurs 
when  the  coordinate  frame  of  the  antenna  (transmitting  or  receiving)  is  different  than  the 
coordinate  frame  of  the  scattering  surface.  Antenna  depolarization  occurs  from  the 
imperfect  isolation  of  the  like-  and  cross-polarized  channels  of  the  antenna.  The  combined 
effect  of  these  two  depolarization  phenomena  can  lead  to  errors  greater  than  the  expected 
cross-polarized  scattering  coefficients  created  by  a  scattering  surface.  Previous  studies  of 
these  effects  have  been  done  in  the  monostatic  case.  This  research  extends  the  previous 
work  to  the  bistatic  case. 

This  paper  is  divided  into  four  separate  sections.  First,  chapter  two  provides  the 
background  of  research  done  in  the  past  to  study  the  scattering  process  and  system 
depolarization  effects.  Chapter  three  provides  the  theoretical  basis  for  extending  known 


system  effects  (translation  depolarization,  antenna  depolarization,  and  scattering 
depolarization)  to  a  bistatic  geometry.  Chapter  four  describes  the  computer  model  designed 
to  simulate  the  bistatic  measurement  of  like-  and  cross-polarization  scattering  coefficients 
from  a  random  surface.  Chapter  five  analyzes  the  results  produced  by  the  computer  model 
to  demonstrate  that  system  induced  depolarization  effects  do  occur  in  the  bistatic 
measurement  environment.  Finally,  chapter  six  presents  the  conclusions  reached  from  the 
analysis  of  simulation  results  and  provides  recommendations  for  decreasing,  or  even 
eliminating,  system  induced  depolarization  effects  from  the  measurement  of  surface 
scattering  in  the  bistatic  environment 


CHAPTER  TWO 


BACKGROUND 

The  complete  use  of  the  polarization  properties  of  scattering  bodies  is  important 
in  order  to  fully  retrieve  all  information  about  a  radar  target.  The  accurate  measurement  of 
all  polarization  characteristics  of  a  target  requires  specific  radar  system  performance  criteria 
in  order  to  assure  correct  data  collection.  This  chapter  presents  the  current  state  of  the 
properties,  descriptors,  and  uses  of  polarization  effects,  as  well  as  the  system  constraints 
necessary  to  retrieve  full  polarization  information. 

In  order  to  fully  understand  the  use  of  polarization  properties  of  radar  targets,  a 
mathematical  basis  must  be  established  to  describe  the  scattering  of  electromagnetic  waves. 
Chan  [1]  reviewed  and  formalized  the  most  common  of  the  polarization  descriptors.  These 
representations  included  time  (and  frequency)  domain,  complex  ratio,  geometric  parameters, 
Stokes  vector,  and  Poincare  sphere.  These  representations  were  extended  to  antenna 
descriptions  and  scattering  matrices.  The  most  general  means  of  identifying  polarization  is 
through  the  use  of  the  polarization  ellipse  with  axes  representing  vertical  and  horizontal 
polarization.  An  example  of  a  polarization  ellipse  is  shown  in  figure  2-1.  All  types  of 
polarization  (i.e.  linear,  circular,  and  elliptical)  can  be  shown  to  be  a  special  case  of 
elliptical  polarization  and  defined  by  a  particular  geometry  on  the  polarization  ellipse.  Four 
important  descriptors  are  needed  to  describe  the  polarization  ellipse,  a  is  the  size  of  the 
ellipse  and  depends  on  the  relative  magnitudes  of  the  wave  in  the  horizontal  and  vertical 
directions.  The  angle  of  tilt,  <)>,  describes  the  tilt  of  the  major  axis  of  the  ellipse  with  respect 
to  an  axis  of  the  coordinate  frame,  x  is  the  ellipticity  of  the  ellipse.  5  is  the  phase  difference, 
in  time,  between  the  orthogonal  components  of  the  coordinate  frame. 
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The  geometric  parameters  (a,  4>,  x,  and  6)  can  be  combined  to  form  Stokes 
parameter  representation  of  polarization.  The  polarization  state  in  Stokes  representation  is 
made  up  of  a  set  of  four  parameters  (g0,  gj,  g2>  g3)  or  (I*  Q»  U,  V}.  Three  of  the  four 
components  are  independent  and  are  related  to  the  geometrical  description  as 


go  =  ka2 

gj  =  k  a2  cos  lx  cos  2<j> 
g2  =  k  a2  cos  lx  sin  2<j) 
g3  =  k  a2  sin  lx 
k  =  (1/2) 1/2 

g02  -  gi2  +  g22  +  S32 


The  Stokes  vector  g  -  [  go»  gj.  g2»  §3  1T  provides  a  unique  descriptor  for  each 
polarization  type  represented  on  the  polarization  ellipse.  For  example,  vertical  polarization 
is  represented  by  g  =  [1,  -1, 0,  0];  horizontal  polarization  by  g  =  [1,1,  0,  0].  Similar 
representations  exist  for  other  linear,  circular,  and  elliptical  polarizations. 

A  graphical  means  of  illustrating  polarization  is  the  Poincare  sphere.  An  example 
of  the  Poincare  sphere  is  shown  in  figure  2-2.  In  the  Poincare  sphere,  normalized  Stokes 
vector  components  are  used  to  map  a  particular  polarization  to  a  unique  point  on  the  sphere. 
Chan  summarized  important  properties  of  use  of  the  Poincare  sphere  to  illustrate 
polarization  [1,  pg.  32]: 

(a)  Each  polarization  ellipse  occupies  a  unique  point  on  the  surface  of  the  sphere. 

The  latitude  and  longitude  of  the  point  are  lx  and  2<j>  respectively. 

(b)  Points  on  the  upper  hemisphere  represent  left-handed  polarization  and 

conversely,  points  on  the  lower  hemisphere  represent  right-handed  polarizations. 

(c)  The  poles  represent  circular  polarizations. 

(d)  All  points  on  the  equator  represent  linear  polarizations,  with  horizontal  (H)  at 

zero  latitude  and  zero  longitude. 

(e)  Orthogonal  polarizations  are  antipodal. 


(f)  All  points  along  a  given  latitude  circle  represent  the  same  ellipticity  (axial 
ratio). 

(g)  All  points  along  a  given  longitude  represent  the  same  orientation  (tilt)  angle 
with  0'  dlt  along  the  zero  meridian  and  90*  tilt  along  the  180°  meridian. 

Chan  showed  the  description  of  scattering  from  an  arbitrary  surface  can  be 
defined  in  terms  of  a  scattering  matrix,  o.  This  matrix  is: 


All  components  of  <T  consist  of  a  magnitude  and  phase.  The  phase  can  be  absolute  or 
relative  to  one  of  the  components.  In  the  bistatic  case,  ct  is  a  function  of  the  transmitter 
incident  angle,  0j,  receiver  incident  angle,  0S,  and  receiver  azimuth  angle,  Os. 

Kennaugh  laid  the  framework  for  the  use  of  polarization  information  in  the  early 
1950’s  with  his  research  in  "Polarization  properties  of  radar  reflections"  [2].  In  the  1970’s, 
Huynen  presented  a  phenomenological  approach  to  the  study  of  polarization  properties  of 
radar  targets  [3-4],  This  showed  the  fundamental  importance  of  connecting  the  geometries 
of  the  target  and  the  transmitting  and  receiving  antennas  in  order  to  fully  describe  the 
electromagnetic  interaction  that  is  taking  place. 

After  the  establishment  of  the  existence  of  important  information  in  polarization 
data,  the  varied  uses  of  this  information  have  been  detailed  by  Boemer  [5-7]  and  Poelman 
[8],  Boemer,  in  particular,  stated  that  high  quality  radar  polarization  data  could  be  used  to 
provide  solutions  to  particular  radar  problems  [8].  These  problems  include: 

(1)  Target  versus  clutter  discrimination:  Polarization  techniques  could 
be  used  to  separate  the  moving  target  from  clutter  existing  in  it's  environment. 
The  clutter  could  be  made  up  of  hostile  chaff  or  surface  effects  close  to  low  flying 
aircraft  The  required  radar  system  would  be  required  to  differentiate  between  the 


relative  dynamic,  unstable,  motion  of  the  target’s  polarization  fork  on  the 
Poincare  sphere  and  the  slowly  changing,  more  stable  polarization  fork  caused  by 
the  clutter. 

(2)  Target-versus-target  and  clutter- versus-clutter  Target  and  clutter 
can  be  separated  by  their  polarization  characteristics  due  to  varying  frequencies. 
Relative  target  and  clutter  size  can  be  discerned  by  increasing  the  broadband 
capabilities  of  the  radar  system  and  observing  the  relative  difference  of  motion  of 
the  polarization  fork. 

(3)  Target  identification:  Complete  collection  of  polarization  and 
doppler  data  from  a  target,  combined  with  electromagnetic  vector  inverse 
processes  and  back-projection  tomographic  techniques  can  be  used  to  uniquely 
identify  and  discriminate  between  targets. 

Research  has  been  done  to  attempt  to  describe  the  actual  process  of  depolarization 
of  electromagnetic  waves  scattered  from  natural  surfaces  [9-15],  The  depolarization 
associated  with  a  radar  target  is  usually  described  as  either  a  surface  scattering  phenomenon, 
a  volume  scattering  phenomenon,  or  a  combination  of  both.  Fung  provided  a  general 
description  of  the  surface  and  volume  scattering  from  a  variety  of  surfaces  [15].  Surface 
scattering  can  be  summarized  [15,  pg.  819]  as: 

Reflection  from  a  smooth-surface  boundary  separating  two  simiinfinite  media  is 
called  specular  reflection  and  is  described  by  the  Fresnel  reflection  laws.  A  wave 
incident  upon  a  rough-surface  boundary  is  partly  reflected  in  the  specular 
direction  and  partly  scattering  in  all  directions.  A  monostatic  radar  (transmitter 
and  receiver  at  the  same  location)  receives  the  backscattered  component  of  the 
scattered  energy.  Thus,  a  monostatic  radar  would,  theoretically,  receive  no 
return  power  from  a  smooth  (specular)  surface  except  for  normal  incidence. 

Therefore,  the  relative  roughness  and  local  incident  angles  determine  how  much 
backscattered  energy  is  returned  to  the  receiver.  The  measurement  of  the  component  of  the 
backscattered  energy  orthogonal  to  the  polarization  of  the  incident  wave  is  needed  to 
characterize  the  depolarization  process  of  the  scattering  surface. 


The  amount  depolarized  scattering  in  a  surface  or  volume  can  be  measured. 
Blanchard  et  al.  [16]  investigated  the  volumetric  depolarization  effects  in  the  returns  of 
cross-polarized  radar  data.  They  concluded  that  target  soil  moisture  was  an  important 
component  in  the  depolarization  process.  However,  they  also  had  to  conclude  that  not  all 
data  analyzed  could  be  used  since  it  was  contaminated  by  the  measurement  process.  Other 
attempts  have  been  made  to  match  theoretical  results  to  experimentally  measured  cross 
polarization  data  [17-19].  Difficulties  have  occurred  when  trying  to  match  the  measured 
cross  section  data  to  the  values  predicted  by  the  theory.  Assuming  that  both  theory  and 
measurements  are  correct  leads  to  the  conclusion  there  must  be  some  phenomena  which 
occurs  during  the  measuring  process  which  contaminates  the  measured  cross  section.  If  the 
phenomenon  can  be  identified  and  applied  to  the  physical  depolarization  measurement,  the 
system  dependency  of  the  measurement  could  be  eliminated. 

One  possible  source  of  error  in  the  measurement  is  antenna  performance  as 
described  by  Blanchard  and  Jean  [22].  They  determined  how  antenna  artifacts  corrupt 
measured  data.  An  important  part  of  antenna  bias  added  to  the  measurement  was  antenna 
feedthrough  of  signals  from  the  orthogonal  input  feed  of  the  antenna.  When  cross 
polarization  information  was  measured,  both  like  and  cross  polarization  channels 
contributed  to  the  total  measured  value,  even  though  only  the  signal  coming  through  the 
cross  polarized  channel  was  desired.  Likewise,  signals  coming  through  the  cross  polarized 
channel  corrupt  the  like  polarization  measurement.  They  described  the  antenna  response 
with  a  matrix  reflecting  the  response  of  each  channel  to  the  orthogonal  wave.  Applying  this 
matrix  to  the  radar  equation  and  solving  for  the  expected  scattering  cross  section, 
demonstrated  the  like- polarized  scattering  coefficient  measurement  showed  very  little 
contamination  from  the  cross- polarized  response  of  the  antenna.  However,  this  was  not  the 
case  for  cross-polarization  scattering  coefficient  measurement.  In  order  to  measure  the 
cross-polarized  scattering  coefficient,  the  receiver's  like-pol  channel  matches  the  cross-pol 


channel  of  the  transmitter.  Feedthrough  of  like-polarized  transmitted  signals  through  the 
cross-polarized  response  of  the  like-polarized  channel  of  the  receiver  significantly  altered  the 
cross-polarized  measurement.  The  isolation  ratio  of  the  like  and  cross  polarized  channels 
was  found  to  be  well  described  at  the  boresight  of  the  antenna.  That  is,  antenna 
manufacturers  commonly  characterize  isolation  ratio  at  boresight,  not  across  the  beamwidth. 
However,  as  the  beamwidth  was  traversed  towards  the  first  null,  the  isolation  ratio  was 
found  to  decrease.  In  a  typical  case,  the  like-  and  cross-pol  responses  could  be  equal  at 
some  point  off-axis.  A  more  proper  way  to  describe  the  isolation  ratio  would  be  to  integrate 
the  isolation  of  the  channels  across  the  entire  beamwidth.  The  isolation  ratio  of  the  antenna 
must  be  large  in  order  to  assure  that  the  measured  scattering  coefficient  is  free  of  any  system 
(antenna)  induced  effects. 

Another  system  induced  effect,  presented  by  Blanchard,  Newton,  and  Jean  [23] 
was  errors  that  could  arise  due  to  the  relative  geometry  of  the  measurement  configuration. 
Electromagnetic  waves  propagated  from  a  transmitting  antenna  are  described,  in  the 
far-field,  as  plane  wave  referenced  to  the  coordinate  system  of  the  antenna.  The  scattering 
surface  is  defined  in  terms  of  it's  own  coordinate  system.  Depolarization  was  found  to 
occur  during  the  transformation  of  the  wave  from  the  antenna  coordinate  system  to  the 
scattering  surface  coordinate  system.  They  found  that  depolarization  was  a  function  of 
beamwidth;  as  the  beamwidth  was  traversed  from  boresight  to  the  first  null,  depolarization 
due  to  translation  increased.  If  the  coordinate  frames  of  the  antenna  and  surface  matched,  as 
can  be  the  case  with  steerable  beam  antennas,  no  depolarization  from  translation  occurred. 
At  nadir,  with  the  antenna  pointing  directly  at  the  surface,  depolarization  occurred  even  at 
boresight.  In  fact,  the  depolarization  effect,  due  to  translation,  at  nadir  was  found  to  be  of 
the  same  magnitude  as  the  cross  polarized  scattering  coefficient  being  measured. 

The  literature  shows  that  research  has  been  done  to  assess  the  effects  of  the 
system  on  measurements.  These  effects  have  been  viewed  primarily  as  to  their  effect  in  a 


monostatic  measurement  environment.  The  effect  of  the  system- induced  depolarization 
effects  in  a  bistatic  measurement  environment  has  not  been  established.  In  a  bistatic 
environment,  less  difference  can  occur  between  cross-  and  like-pol  scattering  coefficients. 
In  order  to  study  the  bistatic  environment,  a  simulation  is  required  which  would  incorporate 
the  known  system-induced  depolarization  effects  into  a  transmitter  and  a  receiver  allowed  to 
move  independently  of  each  other.  It  is  the  purpose  of  this  thesis  to  design  such  a 
simulation  and  extend  the  theory  and  simulation  results  to  the  bistatic  case. 
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CHAPTER  THREE 


ANALYTICAL  DESCRIPTION 


This  chapter  provides  the  theoretical  understanding  and  mathematical  definition 


for  three  depolarization  effects  which  occur  in  a  typical  bistatic  radar  cross  section 


measurement.  The  three  effects  to  be  described  are  translation  depolarization,  antenna 


depolarization,  and  scattering  depolarization.  After  each  effect  has  been  mathematically 


described,  they  will  be  combined  into  a  "system  depolarization  equation"  which  can  be  use 


to  calculate  measured  scattering  coefficients  in  a  computer  model. 


Translation  depolarization  is  the  depolarization  effect  encountered  by  not 


matching  the  coordinate  frame  of  the  transmitting  or  receiving  antenna  to  that  of  the 


scattering  surface.  Translation  depolarization  was  identified  by  Blanchard  et  al.  [23].  In  a 
microwave  measurement  system  all  polarization  vectors  must  be  described  with  respect  to  a 
common,  fixed  reference  point.  The  best,  nonchanging,  reference  point  is  the  scattering 


surface.  Both  receiving  and  transmitting  antennas  are  defined  in  reference  to  the  surface  so 


that  their  relative  interaction  can  be  mathematically  described. 


A  typical  antenna/surface  interaction  geometry  is  shown  in  figure  3-1.  The 
antenna  is  located  in  the  (X'.Y'.Z')  coordinate  frame.  The  surface  is  described  by  the 
(X,Y,Z)  coordinate  frame.  The  origins  of  the  two  coordinate  frames  are  separated  by  a 
distance  h,  defined  by  the  height  of  the  surface  of  the  antenna.  The  Y  and  Y’  axes  are 
always  parallel.  The  antenna  can  rotate  around  the  Y'  axis  in  order  to  point  the  X’  axis  of 


m 


Scattering 


the  antenna  at  a  particular  point.  The  amount  of  rotation  of  the  X'Z’  plane  about  the  Y’  axis 
is  measured  by  the  angle  0O.  80  is  measured  from  the  -Z  axis  to  the  X'  axis.  The  angle  of 
incidence  of  the  individual  incident  wave  is  measured  by  the  angle  0’,  measured  from  the  Z' 
axis.  The  angle  O'  measures  how  far  the  incident  wave  is  off-axis  from  the  main  beam, 
measured  from  the  X'  axis  (O'  =  0°). 

The  antenna  angles,  0'  and  O',  can  be  written  in  terms  of  the  surface  coordinate 
frame  angles,  8  and  O.  In  order  to  solve  for  0,  the  following  intercoordinate  frame 
relationship  is  required: 

0’=  0  +  (y  -  0O) 

thus, 

e  =  e'_  (j-  -  eG) 

Using  this  identity,  0  is  found  from 

cos  0  =  cos  (0  -  (y-80)) 

=  cos  6' sin  0O  +  sin0' cos0o  (3-3) 

sin0  =  sin  (0 ’  —  (y-0o)) 

=  sin0'sin0o  -  cos0' cos0o  (3-4) 

For  the  angle  O,  the  translation  between  two  coplanar  coordinate  frames  is 


(3-1) 

(3-2) 


defined  from  figure  3-2: 


Figure  3-2.  Coplanar  Translation. 


where  (ax,  ay,  *2)  are  unit  vectors  whose  lengths  are  described  in  spherical  coordinates  as: 


|a  x  |  =  sin  0'  cos  G> 
|i  Yl  =  sin  0’  sin  C>' 
I®  Z  I  =  cos0' 


(3-6) 

(3-7) 

(3-8) 


Now,  cos  O  is  defined  in  terms  of  rectangular  coordinates  as: 


a 


X 


Similarly, 


|a  x  |  sin  0  0  +  a  z  cosG  0 

vr|5x|2+  I5  y|2 


sin  9'  cos  C>'  sin60  +  cos  9' cos 9, 

'r|5x|2+  |i  y|2 


sin  <I>  = 


O  xl2  ♦ 


¥  vl 


sin0'  sin  O' 


Ox 


(3-10) 


O  is  calculated  from 


tan  O  = 


sin  9'  sin  O' _ 

sin  0'  cos  <S>  sin0o  +  cos9' cos0o 


(3-11) 


The  magnitude  of  the  incident  wave  can  be  described  as  the  sum  of  two  unit 
vectors,  a'e  and  a'^,  that  propagate  in  the  a’R  direction.  This  vector  must  be  translated  to 
the  surface  coordinated  frame  (unprimed).  The  angle  between  the  incident  wave  and  the 
vertical  plane  of  incidence  on  the  surface  is  defined  as  vf/.  The  relationships  between  4*  and 
the  unit  vectors  (primed  and  unprimed)  in  both  the  0  and  d>  directions  are: 
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O0«O0=O(t)»B<I)  =  COS  ¥ 


(3-12) 


a  e  •  a  (j,  =  -  a  $  •  i  e  =  sin'F 


(3-13) 


These  dot  product  relationships  can  be  solved  with  the  following  spherical  vector 


transformations: 


5$=  -  a  x  sin  O  +  a  Y  cos  O 


(3-14) 


a  <t>  = 


a  'x  sin  Q>  +  5  Y  cos  O' 


=  -(5xsin0o-iz  cos0o)  sin  O  +  a  y  cosO 


=  -  i  x  sin 0O sin  O  +  a  Y  cosO  +  a  z  cos0o  sin  C> 


(3-15) 


The  dot  product  yields 


a  <t»  •  a  d,  =  cos  ¥ 


=  si  i  o  sin  0  sin  d>  +  cos  O  cos  O' 


(3-16) 


A  translation  matrix,  T,  can  now  be  defined  which  completely  describes  the 
transformation  needed  to  get  from  the  antenna  coordinate  frame  to  the  surface  coordinate 


frame.  T  is  defined  as: 


_  P  cos'F  sinT  ”1 

sin  4*  cos  ‘F  J 


(3-17) 


wm. 


The  following  transformation  now  holds: 


i  jjf  '&! 


[::]  •  t[::1 


(3-18) 


The  inverse  relationship,  transforming  from  surface  coordinates  to  the  antenna  frame,  is 
simply  the  inverse  of  equation  (3-18): 


[>:]■ T"  [!:] 


(3-19) 
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Antenna  feed  through  effects,  or  antenna  depolarization,  is  depolarization  caused 
by  imperfect  illumination  of  the  transmitting  or  receiving  antennas.  The  effects  of  antenna 
depolarization  have  been  researched  by  Blanchard  and  Jean  [22].  They  developed  a 
mathematical  basis  for  including  cross  polarization  terms  in  the  radiation  pattern  of  the 
antenna.  Logically,  a  perfect,  or  "ideal",  antenna  would  have  no  coupling  of  vertical  and 
horizontal  transmit  or  receive  modes.  However,  the  real  world  environment  does  not  have 
access  to  perfect  antennas.  Blanchard  and  Jean  characterized  the  term  "polarization  isolation 
ratio"  to  define  the  error  caused  by  cross  polarization  effects  on  a  linear  polarized  antenna. 
They  showed  that  at  boresight  (0°  offaxis)  of  the  antenna  pattern,  the  cross  polarized  return 
had  little  affect  on  the  system.  However,  the  cross  polarized  returns  showed  a  significant 
effect  as  the  pattern  was  traversed  from  boresight  to  the  edge  of  the  pattern.  The 
polarization  response  of  a  typical  antenna  is  shown  in  figure  3-3.  In  this  figure,  the  light 
line  represents  the  antenna  response  for  the  incoming  (or  outgoing)  polarization  of  the  wave 
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i  way  antenna  power  pattern, 
ht  Line  -  Like  polarization  response, 
rk  Line  -  Cross  polarization  response. 


which  matches  the  antenna.  The  dark  line  shows  the  antenna  response  to  the 
orthogonal  wave.  As  can  be  seen,  the  orthogonal  polarization  will  contribute  to  the  off-axis 
response  of  the  antenna.  Blanchard  and  Jean  determined  that  the  polarization  isolation  ratio 
was  a  function  of  the  beamwidth  of  the  antenna.  Cross  polarization  response  added 
considerable  ambiguity  to  measurements  appearing  off  the  main  boresight  axis. 

The  response  of  an  antenna  can  be  described  by  the  matrix  F,  defined  as: 

F  -  [>H  (3-20) 

L1  HV  1 HH  J 

where  the  matrix  terms  are  defined  as  follows: 

fw  =  vertical  response  of  antenna  to  vertical  signal 
fyH  =  vertical  response  of  antenna  to  horizontal  signal 
fHV  =  horizontal  response  of  antenna  to  vertical  signal 
fHH  =  horizontal  response  of  antenna  to  horizontal  signal 

The  matrix  F  represents  the  response  of  the  antenna  to  all  polarizations  for  all  points  in  the 
beamwidth.  In  the  bistatic  environment,  two  F  matrices  are  needed:  FT,  for  the 
transmitting  antenna,  and  F^,  for  the  receiving  antenna. 


Scattering  Depolarization 


A  third  source  of  depolarization  is  the  scattering  surface  itself.  The  scattering 
effect  from  a  surface  seen  through  a  bistatic  geometry  is  significantly  different  than  the 
response  of  the  surface  in  a  monostatic  geometry.  A  scattering  matrix,  o°,  can  be  defined 


0 


to  represent  the  scattering  response  of  the  surface  to  all  polarizations.  o°  is  described  by: 


_  r°w  °vh  "I 
L°hv  °hh  J 


(3-21) 


I 


where 


vertical  scattering  of  vertically  polarized  waves 
vertical  scattering  of  horizontally  polarized  waves 
horizontal  scattering  of  vertically  polarized  waves 
horizontal  scattering  of  horizontally  polarized  waves 


In  order  to  fill  this  matrix,  a  scattering  model  must  be  selected  which  obtains 
practical  results  with  minimum  computational  effort.  Such  a  model  exists  and  is  described 
in  Ulaby  et  al.  [15].  The  model  in  question  is  a  special  case  of  the  Kirchoff  scattering 
model  using  scalar  approximation  and  exponential  correlation.  This  model  exhibits  the 
correct  angular  trend  expected  for  bistatic  scattering  and  provides  results  with  relative 
computational  ease.  The  bistatic  scattering  coefficient,  0°^ ,  is  described  as: 


k2  Kl2 


exp(-q|o2) 


(q2  02)"  -  . 

£  — ^7 -  J  J  p"  exp  [jqx  u  +  j  qy  v]  du  dv  (3-22) 


where 


k  :  wave  number 

a  :  standard  deviation  of  surface  heights 
p  :  surface  correlation  distribution 
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qx  =  k  (sin  0S  cos  0S  -  sin  0T  cos  Ot  ) 
qY  =  k  (sin  0S  sin  0S  -  sin  0T  sin  Ot  ) 
qz  =  k  (cos  0S  +  cos  Ot  ) 


Qz  ~  ^X2  +  ^Y2 


The  surface  correlation  distribution,  p,  is  exponential  in  this  approximation  as  shown 
below: 

p  =  exp  ( -  £  / 1 )  (3-23) 

where 


1  :  correlation  distance  of  surface 
£2  «  u2  +  v2 


Using  this  substitution,  the  integral  in  the  scattering  coefficient  equation  can  be  simplified  to 
a  zeroth  order  Bessel  function.  The  equation  can  then  be  furthered  simplified  to: 


o“  = 
R 


v2  la  I2  ?  ?  00  (q7<J2)nr  _  ■>  ,  ,  -3/2 

21  “P('qZ  °  >  2  [<T )  +  qx  +  qY  ] 


(3-24) 


The  individual  terms  of  the  scattering  matrix  are  found  by  substituting  the  appropriate 
equation  for  aQ  from  the  following  values  of  p  and  q: 


pq  =  HH  : 

a0=_l^  (cos0j+  cos0s  )  cos  (C>s  -d>T ) 

(3-25) 

pq  =  VH  : 

aQ  =~  (  1  +  cos0jcos  0§  )  sin  (<!>$  -  <t>x) 

(3-26) 

pq  =  VV  : 

a0=  r;,  (cos 0x+  cos0 5  )  cos  (<t>s  -Ot  ) 

(3-27) 

pq  =  HV  : 

a0  ="  I^|  (  1  +cos0jcos0s  )  sin  (d>s  -  Oj) 

(3-28) 
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rn  and  Tl  are  the  Fresnel  reflection  coefficients  for  parallel  and  perpendicular 
polarization,  respectively. 


System  Depolarization  Equation 


These  matrices  can  be  combined  in  the  radar  equation  to  assess  their  overall  effect 
on  the  measurement  of  the  scattering  coefficients  of  an  arbitrary  surface.  The  radar  equation 
can  be  rewritten  with  the  matrix  components  as: 


Pr  = 


P  jG jGr X 
(471  )3 


I, 


Ft  T  t 


Ri  Ri 


dA 


(3-29) 


LTOT 


A  system  constant  can  be  defined  to  aid  in  the  normalization  of  this  equation.  The  system 


constant,  Ks,  is 


KS 


PjGj 


gr  x  atot 


(  4  7t )  R  Tnom^  Rnom 


(3-30) 


Ks  and  PR  can  be  combined  to  firm  the  expected,  or  measured,  scattering  coefficient, 
°°M-  This  normalization  results  in: 


(3-31) 


‘TnomKRnom 

A  TOT 


fttt 


Ot  tr 


R 


Ri 


dA 


(3-32) 


Converting  the  integration  over  the  beamwidth  as  described  in  equation  3-32  to  a 
summation  of  the  effect  on  discrete  cells  located  within  the  intersection  of  the  transmitter  and 
receiver  beamwidths  results  in  the  following  equation: 


‘  Tnom1^  Rnom 

A  tot 


2 

cells 


FtTt  ot  tr  fr  aa 

p  2  p  2 

Kt  d 


(3-33) 


The  result  of  this  summation  is  the  approximation  of  the  measured  result  to  the  theoretical 
scattering  coefficients  that  exist  on  the  scattering  surface  with  the  antenna  and  translation 
depolarization  matrices  acting  as  attenuators  or  amplifiers  to  the  surface. 

The  matrix  multiplications  can  be  simplified  to  isolate  the  basic  effect  that 
contribute  to  the  measurement  of  a  like  and  cross  polarization  scatterer.  One  assumption 
that  can  be  made  to  help  this  simplification  is  that: 

cos  *  1  and  sin  *P  *  0  (3-34) 


This  assumption  holds  very  closely  for  all  incidents  angles  except  near  nadir  (0°  incidence). 
With  this  assumption,  the  calculation  of  the  measured  like  and  cross  polarization  scattering 
coefficients  reduces  to  the  following  equations: 


_0 

°HHm  ' 

fRHH 

[aVHfTHV 

+ 

f 

°HH  1 THH 

] 

+  fRHV 

[°VV  f"THV 

+ 

f 

°HVITHH 

] 

(3-35) 
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°HVm  ~ 

fRVV 

[°VHfTHV 

+ 

f 

°HH  ‘THH 

] 

+  fRVH 

[°VV  fTHV 

+ 

CT0  f 
°HV  ‘THH 

] 

(3-36) 

Analysis  of  these  equations  shows  that  the  major  effect  expected  in 
system  depolarization  is  attributable  to  the  cross-polarization  feedthrough  found  in  imperfect 
antennas.  Knowledge  of  the  cross-polarization  feedthrough  of  the  transmitter  and  receiver 
is  necessary  in  order  to  properly  acknowledge  the  source  of  measured  cross  polarized 


scattering  coefficients. 


CHAPTER  FOUR 


MODEL  DESCRIPTION 


A  computer  simulation  was  designed  to  explore  the  combined  effect  of  these 
depolarization  phenomena  associated  with  coordinate  frame  translation,  antenna 
feed-through  effects,  and  scattering.  The  computer  program  was  written  in  FORTRAN  77 
and  was  run  on  a  VAX  11/780.  The  program  simulates  the  surface  interaction  of  a 
transmitter  and  a  receiver,  which  are  oriented  in  an  arbitrary  bistatic  geometry.  The  program 
models  this  interaction  by  taking  a  "snap-shot"  in  time  and  freezing  the  antenna  patterns 
associated  with  both  the  transmitting  and  receiving  antennas.  The  power  received  at  the 
receiver  is  calculated  and  used  to  find  a  measured  value  for  the  scattering  coefficient  of  the 
surface.  This  measured  value  can  then  be  compared  with  the  actual,  or  theoretical,  value  for 
the  scattering  coefficient  which  actually  exists  on  the  surface.  This  chapter  will  describe,  in 
detail,  the  operation  of  the  computer  simulation  as  seen  by  "user". 


is  entered.  The  incremental  spacing  across  the  beam  is  then  input.  The  incremental  spacing 
is  determined  by  the  resolution  of  the  known  data  values  of  the  actual  antenna  pattern  across 
the  beamwidth.  The  user  now  inputs  the  actual  antenna  pattern.  Two  patterns  are  required, 
the  like -polarization  (like-pol)  pattern  and  the  cross-polarization  (cross-pol)  pattern.  The 
user  inputs  the  one-way  power  pattern  values  (in  dB)  found  across  the  pattern  at  9  =  0°  and 
values  of  O  starting  at  0*  up  to  one  half  of  the  beamwidth.  The  incremental  resolution 
determines  the  spacing  of  the  values.  Figure  4-1  illustrates  the  segment  of  the 
antenna  pattern  data  values  input  by  the  user.  The  antenna  pattern  is  assumed  to  be 


Input  Vector 

Lower  quadrant  rotated 
about  input  vector  to  fill 
upper  quadrant 


Figure  4-1.  Antenna  Pattern  Input  Scheme. 

symmetric  about  the  boresight  of  the  antenna  (9  =  0°,  0  =  0°).  This  assumption  is  used  to 
fill  the  lower  quadrant  of  the  pattern  by  rotating  the  input  vector  through  -90°.  Two  value 
linear  interpolation  is  used  to  calculated  values  in  the  the  quadrant  that  are  an  equivalent 
distance  from  the  origin.  The  lower  right  quadrant,  once  it  has  been  completely  filled,  is 


copied  directly  to  the  upper  right  quadrant  to  fill  half  of  the  antenna  pattern.  The  matrix, 
consisting  of  the  upper  and  lower  right  quadrant  of  the  antenna  pattern  for  both  like-  and 
cross-pol,  is  then  written  to  a  file  for  storage. 

Terrain  information  is  also  entered  with  the  DATAINPUT  program.  Two 
different  types  of  terrain  are  used  by  the  model.  The  first  is  monostatic  terrain.  For  this 
type  of  geometry,  the  assumptions  are  made  that  the  terrain  scatterers  act  as  isotropic  point 
sources  and  that  the  receiving  and  transmitting  antennas  are  in  the  same  location.  The 
second  is  bistatic  terrain  Bistatic  terrain  variables  are  assumed  to  change  over  angles  of 
incidence  and  azimuth  in  response  to  the  various  transmitter-to-scatterer-to-receiver 
geometries  that  occur  when  the  transmitter  and  receiver  are  in  different  locations.  Two 
different  subroutines  input  the  scattering  coefficients,  MONOTERRAIN  for  the  monostatic 
environment  and  BISTA'l TERRAIN  for  the  bistatic  environment. 

MONOTERRAIN  allows  the  user  to  enter  the  complex  scattering  coefficient 
matrix,  o°,  for  the  surface.  The  user  enters  a  six  character  code  name  for  the  terrain  file 
about  to  be  entered.  Next,  the  four  complex  components  of  the  a°  matrix  are  entered. 
These  are  ovv,  Oy^,  oHV,  and  o^.  The  numbers  are  entered  as  a  magnitude  (in  dB)  and  a 
phase  (measured  with  respect  to  avv).  Only  one  scattering  matrix  can  be  entered  to 
represent  the  entire  surface. 

The  user  accesses  BISTATTERRAIN  to  input  scattering  coefficient  data  for  a 
bistatic  environment.  The  user  first  enters  a  six  character  code  to  identify  the  data  set  Then 
the  user  enters  the  azimuth  angle  of  the  receiver  and  the  incident  angle  of  the  transmitter. 
The  model  will  only  vary  the  incident  angle  of  the  receiver  for  a  specified  azimuth  and 
transmitter  incident  angle.  Additional  data  files  are  required  in  order  to  store  data  at  other 
receiver  azimuths  or  transmitter  incidents.  The  user  then  enters  a  separate  scattering 
coefficient  matrix,  o°,  for  each  possible  receiver  incident  angle  starting  at  0°,  ending  at  80°, 
incremented  by  10°.  The  four  scattering  coefficient  matrix  components,  as  described  for  the 


MONOTERRAIN  subroutine,  are  entered  as  before. 

The  last  function  of  DATAINPUT  is  to  allow  the  user  to  define  a  "run  set"  for  the 
model  to  operate  on.  A  run  set  is  a  file  containing  all  the  parametric  data  information  needed 
by  the  model  during  it's  execution.  Since  the  model  is  designed  to  be  as  generic  as 
possible,  all  outside  system,  geometry,  and  terrain  information  is  listed  here.  Subroutine 
RUNIN  controls  the  input  of  this  information.  The  user  first  enters  an  eight  character  code 
word  to  identify  this  run  set.  The  user  is  then  asked  two  questions  which  refer  particularly 
to  the  execution  of  the  model.  The  first  is  to  define  matrix  size.  This  size  refers  to  the 
actual  dimension  of  the  matrices  used  determine  discrete  positions  on  the  surface.  The 
higher  dimension  of  the  matrix,  the  greater  the  resolution  of  each  individual  cell  in  the 
matrix,  but  also,  the  greater  execution  time  required  to  run  the  model.  The  next  answer  the 
user  is  required  to  enter  is  the  "X:Y"  ratio.  This  ration  describes  the  relative  size  of  the  cells 
in  the  X  direction  with  respect  to  the  size  of  the  cells  in  the  Y  direction.  Resolution  in  each 
of  the  directions  can  thus  be  changed. 

Antenna  information  is  gathered  next  The  name  of  the  files  containing  the 
antenna  patterns  of  the  transmitter  and  the  receiver  are  entered  by  the  user.  The  user  also 
enters  the  polarization  of  the  antenna,  either  "H"  for  horizontal,  or  "V"  for  vertical. 

The  user  must  now  enter  the  following  information  about  the  relative  positioning 
of  the  receiver  and  transmitter 

(1)  bistatic  or  monostatic  geometry 

(2)  transmitter 

(a)  antenna  inclination  angle,  0O,  in  degrees 

(b)  ground  distance  to  target,  in  meters 

(c)  height  above  surface,  in  meters 

(3)  receiver 

(a)  antenna  inclination  angle 


(b)  slant  distance  to  target 

(c)  minimum  and  maximum  incident  angle 

(d)  azimuth  angle  (transmitter  assumed  to  be  at  CF) 

(4)  terrain  file  name 

After  all  parametric  information  has  been  entered  it  is  stored  in  a  data  file  until  needed  by  the 
execution  of  the  model. 

The  actual  model  is  called  SYSEFF.  This  computer  program  inputs  all 
parametric  information  stored  by  DATAINPUT  and  executes  the  model  accordingly. 
SYSEFF  acts  as  a  stand  alone  program  which  can  execute  without  user  intervention. 

The  first  step  SYSEFF  undertakes  is  to  determine  the  input  file  to  be  used  for 
execution.  The  subroutine  RUNIN  accepts  an  eight  character  code  name  which  is  used  to 
identify  the  input  file.  The  input  file  is  opened  and  all  parametric  variable  values  are  read 
from  this  file.  On  the  basis  of  the  transmitter,  receiver,  and  terrain  file  identifiers  input, 
antenna  patterns  and  scattering  coefficients  are  input  from  their  respective  files.  RUNIN 
then  sets  up  the  looping  mode  of  the  model  to  execute  the  model  the  requested  number  of 
times.  The  architectural  structure  of  the  looping  mode  appears  as: 

(I)  Rx  antenna  inclination  angles 
(A)  Rx  incident  angles 

(1)  Execute  model  (subroutine  RUNONE) 

Subroutine  RUNONE  executes  the  model  one  time  for  a  specific  positioning  of 
the  transmitter  and  receiver.  The  first  operation  of  RUNONE  is  to  "paint",  or  translate,  the 
terrain  with  the  pattern  of  the  transmitting  antenna.  This  is  done  with  the  subroutine 
PAINT.  The  assumption  that  the  antenna  pattern  is  symmetric  through  the  center 
line  defined  by  C>  =  0°  is  used  to  make  these  calculations.  Only  half  of  the  antenna  pattern 
is  "painted"  on  the  ground.  Figure  4-2  shows  a  simple  rendition  of  the  antenna  pattern 
being  translated  to  the  surface. 
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Figure  4-2.  Antenna- to-Surf ace  Pattern  Translation. 

The  surface  is  divided  into  "cells",  each  corresponding  to  a  location  in  the  surface 
matrix.  The  center  of  the  pattern  is  assumed  to  coincide  with  the  center  of  the  surface 
matrix.  PAINT  steps  through  each  cell,  determining  whether  that  particular  cell  is  located 
within  the  beam  of  the  antenna.  If  the  cell  is  located  in  the  beam  of  the  antenna,  the 
following  information  is  calculated  for  that  cell: 

1.  0  and  O  with  respect  to  surface  coordinate  frame 

2.  0'  and  O'  with  respect  to  the  antenna  coordinate  frame 

3.  Like-pol  and  cross-pol  responses  of  the  antenna  (addresses  into 
the  antenna  pattern  matrix  are  0'  and  O',  four  point  interpolation  is  used  to 
calculate  intermediate  values) 

4.  Slant  range  from  the  antenna  to  the  cell 

All  this  information  is  collected  in  a  COMMON  area  for  access  by  other  parts  of  the  model. 

After  the  antenna  pattern  is  "painted"  on  the  ground,  the  pattern  is  folded  along 
it's  access  of  symmetry  in  order  to  complete  the  entire  pattern  on  the  surface.  This  is 


accomplished  with  the  subroutine  DOUBLE. 

Once  the  transmitting  antenna  is  "painted"  and  "doubled"  the  process  continues 
identically  for  the  receiving  antenna.  For  the  receiving  antenna,  however,  an  additional  item 
must  be  taken  under  consideration  following  the  "doubling"  of  the  pattern.  The  transmitter 
is  assumed  to  be  at  0°  azimuth.  The  receiver  can  move  about  in  azimuth  for  a  bistatic 
geometry.  The  receiver  pattern  is  now  rotated  about  the  center  of  the  surface  matrix  to  it's 
correct  azimuth  angle.  This  rotation  is  shown  in  figure  4-3. 


Figure  4-3.  Receiver  Antenna  Rotation. 


After  both  transmitter  antenna  patterns  have  been  placed  on  the  ground  and 
rotated  to  the  correct  orientation,  the  core  of  the  model  is  executed.  The  subroutine 
INTEGRATE  applies  the  radar  equation  derived  in  chapter  three  (equation  3-33)  to  each  and 
every  cell  on  the  surface.  The  five  complex  matrices  (FT,  Tt.  0°,  Tr,  and  FR)  are  all 
known  for  each  cell.  Complex  matrix  multiplication  is  used  to  combine  the  matrices.  Cells 
not  having  both  the  transmitter  and  receiver  pattern  are  skipped.  The  result  of  the 
multiplication,  PR,  is  accumulated  for  each  cell,  resulting  in  a  discrete  integration  of  the 
power  received  over  the  intersection  of  the  transmitter  and  receiver  antenna  patterns.  The 
system  normalization  constant,  Ks,  is  also  calculated  for  each  cell  and  accumulated.  After 
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the  effect  of  all  cells  has  been  included,  the  measured  scattering  coefficient  for  the  surface  is 
calculated  by  dividing  Kg  into  PR. 

Appendix  A  contains  the  following  information  about  DATAINPUT  and 
SYSEFF  useful  to  the  user: 

(1)  computer  listing  of  DATAINPUT 

(2)  computer  listing  of  SYSEFF 

(3)  sample  interactive  session  showing  input  techniques 

(4)  output  results  expected  from  inputs  supplied  at  interactive  session 
This  chapter  described  a  computer  model  designed  to  simulate  the  effect  of 

system  depolarization  effects  in  a  bistatic  environment.  The  model  was  written  in  a  generic 
form  to  allow  the  user  considerable  flexibility  in  the  choice  of  antennas  and  overall 
measuring  system  geometry. 


CHAPTER  FIVE 


RESULTS 

This  chapter  describes  the  experiments  conducted  with  the  SYSEFF  model  and 
the  results  from  these  experiments.  The  experiments  were  designed  specificaHy  to 
determine  the  effect  of  either  translation  or  antenna  depolarization  on  the  bistatic 
measurements  of  scattering  coefficients  on  an  arbitrary  surface. 

The  theoretical  experiments  run  with  SYSEFF,  were  designed  to  determine  a 
measured  scattering  coefficient  value  of  a  specified  surface  as  measured  from  a  variety  of 
bistatic  geometry  positions.  Specifically,  the  positioning  of  the  receiver  and  transmitter 
were  altered  to  discover  what  effect  (if  any)  translation  and  antenna  depolarization  had  on 
measurements  of  like  and  cross  polarized  scattering  from  a  surface.  The  transmitting  and 
receiving  antennas  patterns  were  varied,  as  well,  to  reveal  the  effect  of  the  antenna  patterns 
on  the  measurements. 

The  theoretical  bistatic  surface  scattering  coefficients  used  in  the  actual  test  runs 
were  the  output  of  the  bistatic  scattering  model  developed  by  A.  K.  Fung  [15].  This  model 
was  described  in  chapter  three,  Analytical  Description.  The  bistatic  scattering  coefficients 
were  calculated  for  the  following  bistatic  geometries  in  order  to  establish  the  expected 
scattering  coefficient  matrix  for  each  geometry. 

transmitter  azimuth,  :  0' 

transmitter  incidence,  0T  :  30\60” 

receiver  azimuth,  Os  :  O', 45°, 90°,  1 35\1 80' 

receiver  incidence,  0S  :  0'  -  80' 


normalized  rms  surface  height,  kc  :  0.3 

normalized  surface  height  correlation  length,  kl  :  8 
The  results  of  these  calculations  for  0T  =  60*  are  shown  in  figures  5-1  through  5-3.  The 
results  for  0T  =  30°  are  shown  in  appendix  B.  Each  plot  shows  the  theoretical  scattering 
coefficients  expected  in  the  direction  of  the  scattering  azimuth,  <1>S.  The  following  plots  are 
combined:  0°  and  180°,  45°  and  135°  (225°),  and  ±90°.  A  review  of  the  shapes  of  the 
various  plots  shows  the  following  characteristics  (for  both  30°  and  60°  transmitter  incidence, 
0T).  Strong  forward  scattering  of  like  polarized  components  is  found  on  the  <X>S  =  0°,  180° 
plot  (figure  5-1)  at  approximately  the  angle  as  the  transmitter  incident  angle,  Qj.  No 
scattering  of  cross  polarized  components  is  expected  at  these  scattering  azimuth  angles.  On 
the  d>s  =  45°,  225°  plot  (figure  5-2),  scattering  of  all  four  like  and  cross  polarized 
components  is  predicted.  On  the  Os  =  ±90°  (figure  5-3)  only  cross  polarized  scattering  is 
expected.  This  data  was  used  as  the  scattering  matrices  to  scatter  the  incoming  planes 
waves  from  the  transmitter. 

Three  different  antennas  were  used  for  the  test  runs.  The  transmitter  and  receiver 
were  assumed  to  have  identical  antennas  for  each  run.  The  first  was  an  ideal  antenna, 
having  2.5°  beamwidth  with  0  dB  like  polarization  and  -100  dB  cross  polarization  response 
across  the  beamwidth.  The  second  antenna  was  a  simulated  realization  of  an  actual  antenna. 
The  antenna  pattern  of  this  antenna  is  shown  in  figure  3-3.  The  half-power  beamwidth  of 
this  antenna  was  2.5°.  The  third  antenna  was  a  hybrid  combination  of  the  ideal  and  the  real. 
This  antenna  had  perfect  (0  dB)  like  polarization  response  across  the  2.5°  beamwidth  then 
fell  off  to  constant  -15  dB  sidelobes  out  to  the  full  10°  beamwidth.  The  cross  polarization 
response  of  this  antenna  was  a  constant  -20  dB  across  the  entire  beamwidth.  The  antenna 
patterns  for  the  ideal  and  hybrid  antennas  are  contained  in  appendix  B. 

The  transmitter  was  located  at  a  slant  range  of  10,000  meters  and  an  incident 
angles  of  30°  and  60°  on  the  0°  azimuth.  The  receiver  was  also  located  at  a  slant  range  of 


10,000  meters  but  varied  in  azimuth  and  incidence  angles.  The  receiver  was  moved  through 
0°,  45°,  90°,  135°,  and  180°  azimuth.  At  each  azimuth,  the  receiver  was  taken  through 
0°  -  80°  incident  angles.  Both  transmitter  and  receiver  were  assumed  to  be  pointing  directly 
at  the  same  cell,  thus  as  azimuth  angles  changed,  the  receiver  pattern  on  the  ground  rotated 
around  this  central  cell. 

For  the  measurement  of  the  like  polarization  scattering  coefficient,  both  the 
transmitting  and  receiving  antennas  were  horizontally  polarized.  For  the  cross  polarization 
measurement,  the  transmitting  antenna  was  horizontally  polarized  and  the  receiving  antenna 
was  vertically  polarized.  Both  antennas  were  polarized  with  respect  to  their  respective 
planes  of  incidence. 

The  results  of  the  like  polarization  test  measurements  when  the  transmitter  is  at 
60°  incidence  angle  are  shown  in  figures  5-4  through  5-8.  The  results  for  the  transmitter 
incident  angle  of  30°  are  located  in  appendix  C.  A  comparison  of  the  theoretical  scattering 
coefficients  and  the  measured  scattering  coefficients  for  the  like  polarization  shows  that  little 
significant  difference  can  be  seen  between  the  expected  and  actual  results  except  in  two 
circumstances.  In  figure  5-4  (<DS  =  0°),  all  antennas  measure  the  correct  like-pol  scattering 
coefficient  for  all  angles  except  at  nadir  (0S  =  0°).  At  nadir,  all  antennas  measure 
approximately  3  dB  below  the  expected  scattering  coefficient  Since  all  antennas  show  this 
effect,  it  is  caused  by  the  large  translation  depolarization  that  can  happen  at  nadir  for  even 
very  small  half-power  beamwidths,  in  this  case,  2.5°.  Figure  5-5  (Os  =  45°)  shows  a  very 
similar  curve  to  the  previous  curve.  The  antennas  once  again  were  able  to  measure  the 
like-pol  scattering  coefficient  accurately  for  all  incident  angles  except  nadir.  At  nadir,  the 
ideal  antenna  had  an  error  of  approximately  6  dB.  The  ideal  antenna  had  the  greatest  error 
since  all  of  the  available  power  was  located  within  it's  half-power  beamwidth  and 
depolarized  by  the  translation  effect.  At  90°  azimuth  (figure  5-6),  like-pol  scattering  is 
not  expected.  The  input  for  the  scattering  coefficient  was  set  at  -100  dB.  However,  like 
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Figure  5-6.  Measured  Like-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  90°,  Transmitter  Incidence  =  60°. 


polarized  scattering  was  received.  All  antennas  showed  the  same  range  of  value  (-35  dB  to 
-65  dB)  for  all  receiver  incident  angles.  The  greatest  deviation  of  measured  versus  expected 
like-polarized  returns  occurred  at  C>s  =  135°  (figure  5-7).  This  deviation,  however,  was 
limited  to  between  1  dB  and  2  dB  except  at  nadir,  where  the  effect  of  translation 
depolarization  lowered  all  antenna  responses  to  4  dB  below  the  expected  value.  The 
forward  scatter  plot  (Os  =  180°,  figure  5-8)  shows  perfect  antenna  response  to  the  expected 
like-pol  scattering  response,  except  at  nadir,  where  translation  depolarization  has  it's 
greatest  effect 

Analysis  of  the  cross  polarization  measurements  shows  large  deviations  from  the 
expected  value.  At  both  0“  and  180°  azimuth  (figures  5-9  and  5-10),  no  cross  polarization 
scattering  is  expected,  however  significantly  large  cross  polarization  scattering  was 
measured.  Significant  measured  differences  can  be  seen  directly  attributable  to  antenna 
type.  The  ideal  antenna  came  closest  to  the  expected  cross-pol  values.  It’s  worst 
measurement  was  at  nadir  and  measurements  got  progressively  better  as  incident  angle 
increased.  The  difference  in  expected  and  measured  results  is  caused  by  translation 
depolarization  since  the  ideal  antenna  had  no  cross-pol  feedthrough  to  contribute  to  antenna 
effects.  The  hybrid  and  real  antennas  show  the  added  effects  of  cross  polarization 
feedthrough.  The  real  antenna  at  the  best  measurement  (0S  =  80°)  was  still  60  dB  higher 
than  the  expected  cross-pol  scattering  coefficient.  Close  examination  of  the  forward 
scattering  plot  (<DS  =  180°)  shows  the  measured  cross-pol  return  of  all  antennas  rises  in  a 
manner  similar  to  the  expected  like-pol  scattering  coefficient  at  these  incident  angles  (figure 
5-8).  In  fact,  the  real  antenna  measured  the  cross-pol  coefficient  to  be  almost  exactly  equal 
to  the  theoretical  like-pol  scattering  coefficient  for  all  incident  angles  at  this  azimuth.  This 
indicates  that  the  source  of  the  cross-pol  measurement  is  a  depolarization  phenomena 
inflicted  on  the  like-pol  coefficient.  The  ideal  antenna  shows  the  effects  of  translation 
depolarization.  The  difference  between  the  ideal  and  real  antenna  measurement  curves 
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shows  the  added  effects  of  other  depolarization  phenomena  sources.  Two  sources  are 
immediately  apparent.  The  primary  source  is  antenna  depolarization  from  the  imperfect 
antennas  used  for  transmitting  and  receiving.  Antenna  depolarization  adds  as  much  as  45 
dB  to  the  measurement  at  the  specular  scattering  angle  (  0S  =  60°).  The  secondary  source  is 
additional  beamwidth  that  exists  outside  the  half-power  beamwidth  in  the  real  antenna. 
Translation  depolarization  has  a  greater  effect  the  further  from  boresight  the  antenna  pattern 
occurs.  At  Os  =  45°  and  d>s  =  135°  (figures  5-11  and  5-12,  respectively),  the  ideal  antenna 
measured  the  cross  polarization  scattering  very  accurately  except  at  the  0°  incident  angle, 
where  H*  becomes  large.  The  real  and  hybrid  antennas  showed  significant  differences  in 
measuring  the  cross-polarization  scattering  depending  on  their  relative  cross-polarization 
response.  Once  again  this  difference  is  attributable  to  antenna  feedthrough  of  cross 
polarized  signals.  The  difference  between  the  measured  cross-pol  scattering  coefficient  as 
see  by  the  real  antenna  and  the  ideal  antenna  narrows  as  the  receiver  azimuth  angle 
approaches  90°.  Only  at  the  90°  azimuth  angle  (figure  5-13)  do  all  antennas  measure  the 
cross  polarization  scattering  accurately. 

Significant  results  were  revealed  by  the  theoretical  measurement  the  bistatic 
scattering  coefficients.  These  were: 

(1)  While  measuring  the  like-pol  scattering  coefficient  at  Os  =  90°, 
significant  returns  were  found  where  none  should  have  been. 

(2)  While  measuring  the  cross-pol  scattering  coefficient,  two  important 
results  were  found.  These  were: 

(a)  All  antennas  measured  cross-pol  coefficient  values  much  higher 
than  was  theoretically  expected. 

(b)  The  antennas  themselves  showed  significant  variability  in  the 
measurement  of  the  cross-pol  scattering  coefficient. 


Measured  Cross-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  45°,  Transmitter  Incidence  =  60‘. 
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Receiver  Azimuth  =  90°,  Transmitter  Incidence  =  60”. 


CHAPTER  SIX 


CONCLUSIONS 

This  chapter  discusses  the  conclusions  reached  when  observing  the  results  of  the 
simulation.  Recall,  a  simulation  was  designed  to  model  the  measurement  of  bistatic 
scattering  coefficients.  The  theoretical  scattering  coefficients  used  as  input  to  the  simulation 
were  themselves  theoretical  coefficients  produced  by  a  simplified  Kirchoff  scattering  model. 
Significant  differences  were  found  in  the  comparison  of  the  expected  results  with  the 
simulation  results. 

Two  important  conclusions  can  be  drawn  from  these  results.  One  conclusion  can 
be  drawn  about  translation  depolarization  and  the  second  concerns  antenna  feedthrough. 

First,  the  relative  orientation  of  the  antenna  with  the  ground,  which  results  in 
translation  depolarization,  creates  a  major  depolarization  effect  at  0°  incidence  angle  when 
the  translation  angle,  'F,  becomes  large.  This  effect  is  nullified  when  the  incident  angle  is 
greater  than  10°.  This  effect,  at  nadir,  was  seen  in  all  the  simulation  geometries  tested.  A 
further  contribution  from  translation  depolarization  occurred  when  measuring  cross-pol 
scattering  coefficients  with  the  receiver  azimuth  plane  on  the  transmitter  incident  plane 
(receiver  azimuth  equals  0°  and  180°).  No  cross-pol  scattering  was  expected  in  this 
particular  geometry.  However,  even  the  ideal  antenna,  with  antenna  depolarization 
nullified,  measured  significant  levels  of  cross-pol  scattering. 

The  second  conclusion  that  can  be  drawn  from  the  results,  particularly  the  cross 
polarization  measurements,  is  that  antenna  feedthrough  effects  contribute  a  large  part  to  the 
cross-pol  scattering  measurement  made  at  the  receiver  in  a  bistatic  environment.  This 
conclusion  is  demonstrated  by  the  differences  in  measurement  shown  between  the  real  and 
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ideal  antennas.  The  transmitting  and  receiving  antennas,  suffering  from  cross  polarization 
feedthrough  effects,  deliver  and  accept  extraneous  information  which  can  greatly  increase 
the  measured  value  of  the  cross  polarization  scattering  which  is  expected  off  of  a  given 
surface. 

Clearly,  knowledge  of  translation  and  antenna  depolarization  effects  and  their 
contributions  to  scattering  coefficient  measurements  is  necessary  in  order  to  correctly 
identify  sources  of  cross  polarization  scattering  measurements  in  a  bistatic  environment. 
Three  areas  of  the  measurement  system  must  be  controlled  in  order  to  minimize 
measurement  of  system  induced  depolarization  and  maximize  measurement  of  target  induced 
depolarization.  These  are: 

(1)  No  measurements  should  be  taken  at  or  near  nadir.  Significant 
translation  depolarization  occurs  at  nadir  even  when  using  "ideal"  antennas. 

(2)  Antennas  should  be  chosen  to  assure  the  maximum  isolation  ratio 
between  like  and  cross  polarized  channels  across  the  entire  beamwidth  of  the 
antenna. 

(3)  Antennas  should  be  chosen  with  minimum  beamwidth  in  order  to 


nullify  translation  depolarization  effects  which  occur  off-axis. 
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0001 

PROGRAM  DATA1.IPOT 

0002 

- 

0003 

2*  *  * 

•  •• 

*  * 

0004 

2 

THE  PURPOSE  Or  T.U  -  PROGRAM  IS  TO  INPUT  THE 

DATA  NEEDED  TO 

0005 

z 

CALCULATE  THE  VALUES  OF  RADAR  BACX3CATTERISG 

COEFFICIENTS. 

0006 

2 

DATA  CONCERN l 'IG  ANTENNA  LIKE-  AND  CROSS-POLARIZATION  VALUES, 

0007 

z 

AND  TERRAIN  CHARACTERISTICS  IS  PLACED  INTO  APPROPRIATE  FILES. 

0303 

2 

the  program  is  then  used  to  make  the  desired 

CALCULATIONS, 

0009 

w 

USING  THE  DAT-.  INPUT  HERE. 

• 

0010 

2  +  *  + 

k  * 

0011 

z 

0012 

2 

SET  UP  SELECTION  MENU  FOR  USER 

0013 

2 

0014 

10 

WRITE (6 ,700 ) 

0015 

700 

FORMAT  Cl’ ,////, 25X,* SYSTEM  EFFECTS  ANALYSIS 

• 

9 

0016 

1 

///,’  1.  INPUT  ANTENNA  DATA’, 

0017 

2 

/,'  2.  INPUT  RUN  SET  DATA’, 

0018 

3 

/,’  3.  INPUT  I30TR0PIC/M0N0STATIC  TERRAIN 

FILE’  , 

0019 

4 

/,’  4.  INPUT  BISTATIC  TERRAIN  FILE’, 

0020 

5 

/,'  5.  END1) 

0021 

n 

0022 

ICHOICE-GETVAR  ( ’  ENTER  CHOICE  -  •) 

0023 

c 

0024 

2 

FOR  AN  INCORRECT  CHOICE,  RETURN  TO  SELECTION 

MENU. 

0025 

c 

0025 

IF  ( ICHOICE.LT . 1 .OR . ICHOICE.GT . 5)  GOTO  10 

0027 

GO  TO  (50,100,150,200,999)  ICHOICE 

0028 

2 

0029 

50 

CALL  ANTIN 

0030 

GO  TO  10 

0031 

2 

0032 

100 

CALL  RUN  IN 

0033 

GO  TO  10 

0034 

C 

0035 

150 

CALL  MONOTERRAIN 

0036 

3D  TO  10 

0037 

C 

0033 

200 

CALL  3ISTATTERRAIN 

0039 

GO  TO  10 

0040 

c 

0041 

999 

CONTINUE 

0042 

2 

0043 

END 
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3001 

SUBROUTINE  ANT IN 

'->002 

2 

0003 

*  * 

*** 

0004 

2 

INPUT  ANTENNA  PATTERN  DATA 

0005 

^  •  * 

*  *  • 

000c 

z 

0007 

DIMENSION  LIKE (-50:50,0:50) , CROSS (-50:50,0:50) 

0008 

REAL  LIKE 

0009 

CHARACTER'S  ANT 

0010 

CHARACTER* 1 3  ANTFILNAM 

0011 

LOGICAL  AGAIN, YESNO 

0012 

COMMON  /POINTS/BWINC.MAXPT 

0013 

2 

0014 

WRITE (6,700) 

0015 

700 

FORMAT (' 1 ’,////, 25X, ’ANTENNA  DATABASE  INPUT  ROUTINE*) 

0016 

2 

0017 

c 

CALL  SUBROUTINE  ’GETANTNAM'  TO  INPUT  ANTENNA  FILE  NAME. 

0018 

2 

IF  FILE  ALREADY  EXISTS,  ASK  IF  IT  IS  TO  BE  REPLACED. 

0019 

2 

0020 

CALL  GETANTNAM (ANT .ANTFILNAM, STATUS) 

0021 

IF  (STATUS .EO . 1 )  THEN 

0022 

AGAIN-YESNOC  REPLACE  (Y  ,N)  ?’ ) 

0023 

IF  ( .NOT. AGAIN )  GO  TO  999 

0024 

ENDIF 

0025 

z 

0026 

2 

ENTER  BEAMWIDTH  AND  INCREMENTAL  RESOLUTION  ACROSa  BEAMvJIDTH 

0027 

z 

VALUES  ARE  ENTERED  IN  DEGREES 

0028 

n 

0029 

BW-GETVARC  ENTER  3EAMWIDTH :' ) 

0030 

BWINC-GETVAKC  ENTER  INCREMENTAL  RESOLUTION  ACROSS  BEAM'WIDTH:  '  ) 

0031 

c 

0032 

c 

DEFINE  NUMBER  OF  DATA  POINTS  ACROSS  B3AMWIDTH 

0033 

c 

0034 

BWPTS-BW/BWINC 

0035 

IF ( (MOD (BWPTS ,2.0) ) .EQ.O)  BWPTS-BWPTS+i 

0036 

MAXPT-INT (BWPTS/2 ) +1 

0037 

c 

0038 

WRITE (6,702) 

0039 

702 

FORMAT (/,'  •*•  ENTER  LIKE  POLARIZATION  RESPONSE’, 

0040 

1 

'  ACROSS  BEAMWIDTH  ***’,//) 

0041 

CALL  GETARRAY (LIKE) 

0042 

CALL  FILL (LIKE) 

0043 

C 

0044 

WRITE (6,704) 

0045 

704 

FORMAT (/,'  •*•  ENTER  CROSS  POLARIZATION  RESPONSE’, 

0046 

1 

'  ACROSS  BEAMWIDTH  ***',//) 

0047 

0048 

CALL  GETARRAY (CROSS) 

0049 

CALL  FILL (CROSS) 

0050 

c 

0051 

c 

PASS  DATA  TO  SUBROUTINE  FOR  STORAGE 

0052 

2 

0053 

CALL  STORANT (ANT , ANTFILNAM , BW , BWINC , MAXPT , L I KE .CROSS ) 

0054 

n 

0055 

99? 

CONTINUE 

0056 

c 

005-> 

RETURN 

0058 

END 
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0001 

o..'  >2 

0003 

0004 

0005 

0006 

0007 

0003 

0009 

uiK 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

002b 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

005" 

0058 

0059 

0060 


SU3.40UTIN2  RUN IN 

A  ft***#*******************************.* ~  .<  *  • 

INPUT  RUN  SET  DATA  * 


COMMON  /LARGE/S l  ZS , YS I ZS , ARE A 
CHARACTER*!*  TXA.4  T  ,  RXANT  ,  TERN  AM 
CHARACTER*8  RUNNAM 
CHARACTER *11  RUNFILNAM 

CHARACTER* 10  TXFILNAM , RXFILMAM .TSRFILNAM 
CHARACTER *1  MONOSTATIC , RXNOMOMATCrf ,ANTPOLAR 
INTEGER  STATUS , FX VSRPOL , TXHORPOL , RXVERPOL , RXHORPOL 

C  INPUT  NAME  OF  RUN  FILE.  CREATE  AN  INPUT  FILE  WITH  GIVEN  NAME 

WRITS (6,800) 

300  FORMAT (' 1 1 ,////, 25X , *  RUN  SET  DATA  INPUT  ROUTINE’) 

WRITS (6,700) 

700  FORMAT (//,’  INPUT  NAME  OF  RUN  FILS') 

ACCEPT  701, RUNNAM 

701  FORMAT (A8) 

RUNFILNAM-RUNNAM//’ .IN’ 

SIZS-GETVARC  ENTER  MATRIX  SIZE’) 

YSIZE-G3TVARC  ENTER  X:Y  RATIO  (1:?)') 

C  GET  T  <  ANTENNA  NAME 

WRITE (6,702) 

702  FORMAT (/,’  •**  TRANSMITTER  ***',/) 

CALL  GETANTNAM (TXANT , TXFILNAM , STATUS ) 

IF  (STATUS. EQ.O)  THEN  !  ABORT  PROGRAM  IF  FILE  DOESN’T  EXIST 
WRITE (6,703) 

703  FORMAT (’  RUN  ABORTED  **•  TX  ANTENNA  FILE  NOT  FOUND') 
GOTO  999 

ENDIF 

C  GET  TX  ANTENNA  POLARIZATION 

205  WRITE (6 , 300 ) 

300  FORMAT (/ , '  ENTER  TRANSMITTER  POLARIZATION  (V,H)-') 

n 

ACCEPT  707  ,AHTPOLAR 


IF  (ANTPOLAR . EQ . ’ V ’ )  THEN 
TXVERPOL-1 
TXHORPOL=0 

ELSE 

IF  (ANTPOLAR. EQ. ’H*)  THEN 
TXVERPOL=0 
TXHORPOL*=l 

ELSE 

GOTO  205 

ENDIF 

ENDIF 
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RUN  IN 
0061 

- 

0062 

GET  RX  ANTENNA  NAME 

0063 

2 

0064 

WRITE  (6 ,704 ; 

0065 

704 

FORMAT (/,’  ***  RECEIVER  ***',/) 

0065 

z 

0067 

CALL  G2TANTN AM ( RXAX  T  ,  RXF I LN  AM , STATUS ) 

006? 

I?  15TATU3 .  SJ.O)  THEN  ;  ABORT  RUN  IF  FILS  DOES..’ '  T  i 

0069 

WRITE (6,705) 

007.) 

"05 

FORMAT ( '  RUM  ABORTED  ***  RX  ANTENNA  FILE  NOT  FOUND’) 

0071 

GOTO  999 

0072 

END  IF 

0073 

Z 

0074 

2 

GET  RECEIVER  ANTENNA  POLARIZATION 

0075 

2 

0076 

210 

WRITE (6,310) 

0077 

310 

FORMAT (/,'  ENTER  RECEIVER  POLARIZATION  (V , H ) -  * ) 

0073 

ACCEPT  707 , ANTPOLAK 

0074 

0080 

IF  (ANTPOLAR.EQ. 'V' )  THEN 

0031 

RXVERPOL-1 

0082 

RXHORPOL-O 

0083 

ELSE 

0084 

IF  (ANTPOLAR.EQ. 'H' )  THEN 

0085 

RXVERPOL-O 

0086 

RXHORPOL-1 

0037 

ELSE 

0088 

GOTO  210 

0089 

ENDIF 

0090 

ENDIF 

0091 

w 

0092 

WRITE (6,706) 

0093 

706 

FORMAT (/,'  RECEIVER  AND  TRANSMITTER  IN  SAME  LOCATt')*  (Y,N)?’) 

0094 

W 

0095 

ACCEPT  707, MONOSTATIC 

0096 

707 

FORMAT (Al) 

0097 

C 

0093 

2 

GET  TRANSMITTER  POSITION  IF  TRANSMITTER,  RECEIVE1*  AKIJ  NOT  IN 

0099 

z 

THE  SAME  LOCATION 

0100 

c 

0101 

IF  (MONOSTATIC. NE. *Y’ .AND. MONOSTATIC. NE. 'y' )  THEN 

0102 

TXMID-GETVAR ( '  ENTER  GROUND  DISTANCE  FROM  TX  TO  TARGET ' 

0103 

TXHEIGHT-GETVARC  ENTER  TRANSMITTER  HEIGHT’) 

0104 

ENDIF 

0105 

c 

0106 

RXSLANT-GETVAR (’  ENTER  SLANT  DISTANCE  FROM  RX  TO  TARGET ') 

0107 

RXNOMTHETAMIN-GETVAR ( '  ENTER  MINIMUM  RX  INCIDENT  ANGLE') 

0108 

RXNOMTHETAMAX-GETVAK ( '  ENTER  MAXIMUM  RX  INCIDENT  ANGLE') 

0109 

c 

0110 

c 

ASK  IF  RECEIVER  ANTENNA  INCLINATION  MATCHES  INCIDENT  ANGLE. 

Oil  L 

z 

IF  NOT,  GET  MINIMUM  AND  MAXIMUM  RECEIVER  ANTENNA  INCLINATIONS. 

0112 

z 

0113 

WRITE (6 ,708 ) 

0114 

708 

FORMAT (/, '  RX  ANTENNA  INCLINATION  MATCH  INCIDENT  ANGLE  (Y ,N  )  ?  ' 

0115 

z 

0116 

ACCEPT  709 , RXNOMOMATCH 

0117 

709 

FORMAT (Al) 

59 


run:-. 

.  -i 

0119 

IF  <RX!JOM9MATCH.NE. *Y' .AMO . RXNOMO MATCH .SE . ’v’  I  THEM 

0120 

RXTH2TA0MIN-GSTVAR { '  ENTER  MINIMUM  RX  ASTSSSA  INCLINATION 

0121 

RXTHSTA0MA\=32rVAR  ( '  ENTER  MAXIMUM  RX  ANTENNA  ISCLI SAT I UN 

.  •  *  0 

ELSE 

0123 

RXTHETAOMIN-RXNOMTHETAMIM 

0124 

RXTHETAOMAX-RXNOMTHETA.MA  . 

0125 

ENDIF 

0126 

I— 

012? 

RXAZI M-GETVAR ( '  ENTER  RX  AZIMUTH  ANGLE’) 

0123 

Sm 

0129 

2 

GET  DESIRED  TERRAIN  FILE.  ABORT  RUN  IF  FILE  DOESN’T  EXIST 

013U 

£ 

0131 

CALL  GET TERN AM ( IHRSAM ,TERFI LNAM ,STATUS ) 

0132 

IF  (STATUS. EQ.O)  THEN 

0133 

WRITE  (6,715) 

0134 

715 

FORMAT ('  RUN  ABORTED  TERRAIN  FILE  NOT  FOUND') 

0135 

GOTO  999 

0136 

ENDIF 

0137 

c 

0133 

c 

OPES  RUN  FILE  AND  ENTER  INPUT  DATA 

0134 

0140 

OPEN (UNIT-20, FILE-RUNFILNAM, STATUS- ’NEW' , RECORDS  I ZE-1000 ) 

0141 

2 

0142 

WRITE ( 20 , 200 ) MONOSTATIC , RXNOMOMATCH , TXANT , RXANT , 

0143 

1 

TXMID ,TXHE IGHT , 

0144 

2 

RXSLANT , RXMOMTHETAMI H , RXNOMTHETAMAX , 

0145 

3 

RXTHETAOMIN , RXTHETAOMAX , 

0146 

4 

RXAZtM,SIZE,YSIZE,TERNAM, 

0147 

5 

TXVERPOL , TXHORPOL , RXVERPOL , RXHORPOL 

0143 

200 

FORMAT <2A1,2A6,10F14. 7, A6, 411) 

0149 

c 

0150 

CLOSE (UNIT-20) 

0151 

c 

0152 

c 

CREATE  COMMAND  FILS  WITH  RUN  FILS  NAME  TO  RUN  PROGRAM  'SYSEFF 

0153 

c 

0154 

OPEN  (UNIT-20, FILE-RUNNAM//' .COM' , STATUS- ’ NEW’ ) 

0155 

WRITE (20,201) RUNNAM 

0156 

201 

FORMAT  (’$  RUN  SYSEFF’ ,/,A8) 

0157 

CLOSE  (UNIT-20) 

0158 

j 

0159 

999 

RETURN 

0160 

END 
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C  v  v.  _ 
0002 

0003 
0  j04 
0005 
0006 
0007 
000  5 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0013 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 


SUBROUTINE  MONOTERRAIN 


INPUT  MONOS T A. TIE  TERRAIN  3ATA 


CHARACTSR*6  TER 
CHARACTER* 10  TERFILNAM 
LOGICAL  AGAIN , YSSNO 

WRITE (6,700} 

700  FORMAT (' 1 ',////, 25X, 

1  ' MONOSTATIC  TERRAIN  D  AT  Si  L  IN?  ■'  :  C  *  ) 

C  GET  A  SIX-CHARACTER  TERRAIN  FIuE  IDENTIFIER.  ASK  IF  ANY 

C  EXISTING  VERSION  OF  Ttii  FILE  SHOULD  BY  REPLACED. 

CALL  GETTERNAM (TER , TERFILNAM , STATUS ) 

IF  (STATUS. EQ.l)  THEN 

AGAIN-YE3N0C  REPLACE (Y , N )?' ) 

IF  (.NOT. AGAIN)  GO  TO  999 

ENDIF 

C 

C  GET  SCATTERING  COEFFICIENTS 

SIGW-GETVARC  ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  '  ) 
ANGW-GETVAR  ( '  ENTER  V7  PHASE  ANGLE  '  ) 

SIGVH-GETVARC  ENTER  VH  SCATTERING  COEFFICIENT  UN  DB) ’ ) 
ANGVH-GETVAR ( '  ENTER  VH  PHASE  ANGLE') 

SIGHV-GETVAR ( '  ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)') 
ANGHV-GETVAR ( '  ENTER  HV  PHASE  ANGLE’) 

5IGHH«GETVAR ( '  ENTER  HH  SCATTERING  COEFFICIENT  (IN  03)') 
ANGHH-GETVAR  ( '  ENTER  HH  PHASE  ANGLE') 

C  OPEN  DATA  STORAGE  FILE 

OPEN  (UNIT»20 ,  FILE «TERF I LN AM .STATUS* ' NEW ' .REC0RD5I ZE-500) 

C 

WRITE  (20,200) SI GW ,SIGVH , S IGHV , SIGHH , ANGVV  ,  ANGVH , ANGHV , ANGHH 
200  FORMAT  (8F14.7) 

C 

CLOSE  (UNIT=20) 

n 

999  CONTINUE 

C 

RETURN 

END 
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3001 

0002 

2  *  •  * ' 

0003 

z 

3304 

2**  * 1 

0003 

0006 

0007 

0308 

0009 

0010 

z 

0OU 

0012 

700 

0013 

2 

0014 

z 

0015 

z 

0016 

2 

0017 

0013 

0019 

0020 

0021 

0022 

z 

0023 

z 

0024 

2 

0025 

0026 

0027 

2 

0028 

2 

0029 

2 

0030 

z 

0031 

0032 

0033 

0034 

701 

0035 

z 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

50 

0046 

n 

0047 

z 

0048 

z 

0049 

0050 

0051 

199 

0052 

z 

0053 

0054 

60 

0055 

0056 

200 

0057 

0058 

999 

0059 

0060 

SUBROUTINE  3 1  ST  ATT  E  RR  AI  .4 


INPUT  3 1. STATIC  TERRA  I  :  DA;’-. 


DIMENSION  S IGVV (0:8) ,5IGVH(0:3) ,3IGHV(0:8) ,SIGHH(0:3) 

DIMENSION  ANGV7 (0:3) ,ANGVH(0:3) ,ANGH7(3:3) ,ANGHH(0:3) 
CHARACTER'S  TER 
CHARACTER *10  TERFILNAM 
LOGICAL  AGAIN, YESNO 

i  JR ITE(6,700) 

FORMAT ( ’ 1  * ,//// ,  25X , ’ BISTATIC  TERRAIN  DATA  FILE  INPUT  ROUTINE’) 

GET  A  SIX-CHARACTER  TERRAIN  FILE  IDENTIFIER.  ASK  IF  AN Y 
EXISTING  VERSION  OF  THE  FILE  SHOULD  BE  REPLACED. 

CALL  GETTERNAM (TER , TERF I LNAM , STATUS ) 

IF  (STATUS. EQ.l)  THEN 

AGAIN-YE3N0C  REPLACE  (X  ,  N )?  1  ) 

IF  (.NOT. AGAIN)  GO  TO  999 

END  IF 

ENTER  RECEIVER  AZIMUTH  ANGLE,  TRANSMITTER  INCIDENT  ANGLE 

AZIM-GETVAR ( 1  ENTER  RX  AZIMUTH  ANGLE’) 

TXINC-GETVAR ( ’  ENTER  TX  INCIDENT  ANGLE’) 

FOR  INCIDENT  ANGLES  PHI  »  0  TO  80  DEGREES,  GET  SCATTERING 
COEFFICIENTS  AND  PHASE  ANGLES 


DO  50  I >0 , 8 
?HI-10*I 

WRITE (6 ,701 )  PHI 

FORMAT ( ’  INCIDENT  ANGLE-’ ,F5. 2) 


3IGW  (I )  -GETVAR  ( ’ 
ANGW  ( I )  -GETVAR  ( ’ 
SIGVH (I ) -GETVAR  (’ 
ANGVH (I ) -GETVAR ( ’ 
SIGHV (I ) -GETVAR (’ 
ANGHV (I ) -GETVAR (’ 
3IGHH (I ) -GETVAR  (’ 
ANGHH  (I ) -GETVAR  (’ 


ENTER  W 
ENTER  W 
ENTER  VH 
ENTER  VH 
ENTER  HV 
ENTER  HV 
ENTER  HH 
ENTER  HH 


SCATTERING  COEFFICIENT  (IN  DB)’) 
PHASE  ANGLE’) 

SCATTERING  COEFFICIENT  (IN  DB)  ’  ) 
PHASE  ANGLE’) 

SCATTERING  COEFFICIENT  (IN  DB)  ’  ) 
PHASE  ANGLE’) 

SCATTERING  COEFFICIENT  (IN  DB)  ’  ) 
PHASE  ANGLE’) 


CONTINUE 


OPEN  TERRAIN  DATA  STORAGE  FILE 


OPEN  (UNI T-20 , FILE-TERFI LNAM , STATUS- ’ NEW ’ , RECORDS  I ZE-S00 ) 
WRITE  (20,199)  AZIM.TXINC 
FORMAT  (2F8.3) 

DO  60  1-0,8 

WRITE  (20,200) SIGW  ( I )  ,SI3VH (I )  , SIGHV (I )  , S IGHH ( I )  , 

ANGW  (I  )  , ANGVH  (I )  , ANGHV  (I )  , ANGHH  (I  ) 

FORMAT  (8F14.7) 

CLOSE  (UNIT-20) 

CONTINUE 

RETURN 

END 


3.X  . 

3032 

OOOi 

z 

SUBROUTINE  GETAr.-  i^N.A.M  (ANT , AN FFILNAM, STATUS ) 

-  *  *  *  " 

0004 

2 

input  antenna  identifier 

»Vi  *5 

0006 

0007 

CHARACTER*6  ANT 

0003 

CHARACTER*!..  ANTFILNAM 

0009 

z 

0010 

2 

GET  SIX-CHARACTER  NAME  FOR  ANTENNA  FILE 

0011 

2 

0012 

WRITE (6,700} 

Ool  3 

700 

FORMAT (’  ENTER  ANTENNA  IDENTIFIER:') 

00t4 

0015 

ACCEPT  701, ANT 

0016 

701 

FORMAT (A6) 

0017 

2 

0018 

ANTFILNAM-ANT// ' .ANT ' 

0019 

2 

0020 

2 

OPEN  FILE  WITH  GIVEN  NAME;  INDICATE  IF  FILT  ALR'.VK 

0021 

c 

0022 

OPEN (UNIT-10,  FILE-ANTFIL.MAM, STATUS-' OLD' , ERR=10 ) 

0023 

CLOSE (UNIT-10) 

0024 

z 

0025 

WRITE (6,702) 

0026 

702 

FORMAT {’  FILE  FOUND') 

0027 

STATUS-1 

0028 

GOTO  99 

0029 

2 

0030 

10 

WRITE (6,703) 

0031 

703 

FORMAT ('  FILE  NOT  FOUND') 

0032 

0033 

STATUS-0 

0034 

0035 

99 

RETURN 

0036 

END 
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0001 

subroutine  ge ttsrnam ctsr, terfiln am, status ) 

0002 

z 

0002 

z  *  *  9  * 1 

0004 

z 

1.4?  TERRAIN  IDENTIFIER 

0005 

•  •  •  1 

0006 

c 

0007 

CHARACTER* 5  TE- 

■  i008 

CHARACTER*  10  TSRFILNAM 

0009 

2 

0010 

z 

GE:’  5 i {-CHARACTER  TERRAIN  FILS  NAM? 

0011 

2 

0012 

WRITE(6,700) 

0013 

700 

FORMAT  ENTER  TERRAIN  FILS  IDENTIFIER:’) 

0014 

c 

0015 

ACCEPT  701, TER 

0016 

701 

FORMAT (Ao) 

0017 

z 

0018 

TERFILNAM-TER//' .TER' 

0019 

2 

0020 

2 

OPEN  FILE  WITH  GIVEN  NAME;  INDICATE  IF  FILS  ALREADY  EXISTS 

0021 

2 

0022 

OPEN (UN IT* 10 , F ILE -TERF I LNAM , STATUS- ' OLD ' ,ERR=10) 

0023 

CLOSE (UNIT-10) 

0024 

2 

0025 

WRITS (6,702) 

0026 

702 

FORMAT (’  FILE  FOUND’) 

0027 

STATUS-1 

0023 

GOTO  99 

0029 

C 

0030 

If. 

WRITE (6, 70 3) 

0031 

703 

FORMAT (’  FILE  NOT  FOUND’) 

0032 

STATUS-0 

0033 

r* 

0034 

99 

RETURN 

0035 

END 
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0031 

SUBROUTINE  GETARKAY (MAT) 

0002 

2 

00  * 

3004 

2 

INPUT  ARBITRARY  ARRAY  OF  ANTENNA  PATTERN  DATA 

# 

0005 

00  ou 

2 

3007 

REAL  MAT (-50: 50, 0:50; 

000d 

COMMON  /POIMTS/BWINC , MAXPT 

0009 

2 

COLO 

2 

INPUT  ONE-WAY  POWER  PATTERN  VALUES  (IN  DB)  FOR  EACH 

0011 

2 

INCREMENTAL  ANGLE 

0012 

2 

0013 

DO  10  I “0 , MAXPT— 1 

0014 

PHI»HWt:lC»I 

0015 

WRI TS (6 , 700 )  PHI 

- 

0016 

700 

FORMAT c  ANGLE’ ■  ,F5.2,  1  DB  = ' ) 

0017 

READ (5 , * )  MAT (0,1) 

0018 

2 

0019 

2 

CONVERT  FROM  DB 

0020 

0021 

MAT (0,1) »10«* (MAT (0,1) /20 ) 

0022 

r* 

V. 

0023 

lf. 

CONTINUE 

0024 

2 

0025 

RETURN 

0026 

END 

0001 

SUBROUTINE  FILL (MAT) 

0002 

2 

0003 

£*##* 

0004 

c 

FILL  ANTENNA  PATTERN  MATRIX  WITH  VALJES  VkSED  UPON  INPUT  VECTOR 

0005 

0006 

w 

0007 

REAL  MAT (-50:50,0:50) ,MAXDIS ,PNTDIS 

0008 

COMMON  /POINT3/BWINC, MAXPT 

0009 

2 

0010 

MAXDIS’DIST (0 .MAXPT-l) 

0011 

2 

0012 

DO  10  I«l, MAXPT-l 

0013 

DO  10  J-0, MAXPT-l 

0014 

PNTDIS-DIST (I ,J) 

0015 

C 

0016 

2 

INTERPOLATE  BETWEEN  KNOWN  INPUT  POINTS 

0017 

0018 

IF(PNTDIS.GT.MAXDIS)  THEN 

0019 

MAT  (I ,  J  )  -  0 

0020 

ELSE 

0021 

II-ABS (INT (PNTDIS) ) 

0022 

FRAC-PNTDIS-II 

0023 

DIFF-MAT (0,11+1) -MAT (0,11) 

0024 

MAT ( I , J ) -MAT (0,11) +FRAC  *D I FF 

0025 

ENDIF 

0026 

2 

0027 

MAT (-1 ,J) -MAT (I ,J) 

0028 

c 

0029 

i) 

CONTINUE 

0030 

*2 

0031 

RETURN 

0032 

END 

0001 

SUBROUTINE  GET ANT (ANT , NAM ,3W, BWINC , MAXPT 

0003 

j  j04 

c 

READS  ANTENNA  OAT  .-’RDM  FILE 

0005 

2  *  *  *  * 1 

j0  06 

2 

0007 

DIMENSION  LIKE (-50:50,0:50)  , CROSS (-50 : 50 

0008 

REAL  LIKE 

0007 

CHARACTER*6  ANT 

0010 

C H A R AC T E  P.  *  1 0  NAM 

0011 

z 

0012 

2 

OPEN  A  FILE  WITH  ANTENNA  IDENTIFIER  NAME 

0013 

2 

0N1 

>PSN  (UNIT-10 , F ILE-NAM , STATUS* 1  OLD ' ) 

0015 

2 

0016 

READ (10,700)  ANT, BW, BWINC, MAXPT 

0017 

■;o 

FORMAT (A6 , 2F3 .3,14) 

0013 

2 

0019 

DO  10  I«- (MAXPT-1) , MAXPT-1 

0020 

10 

READ (10,701) (LIKE ( I , J ) , J-0 , MAXPT-1 ) 

0021 

DO  20  I  a- (MAXPT-1 ), MAXPT-1 

0022 

20 

READ (10, 7  1) (CROSS ( I, J) , J-0 , MAXPT-1 ) 

0023 

701 

FORMAT (100F3. 3) 

0024 

2 

0ii25 

CLOSE (UNIT-10) 

0026 

2 

0027 

RETURN 

0023 

END 

0  0  0  L 

0002 

0003 

0004 

0005 

0006 

0007 

0003 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0013 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 


.5'JBROUTINE  3TORAMT  (ANT  , NAM, BW, BWINC  ,MAXPT  , LI  KE , CROSS ) 


STORES  ANTENNA  IJATA  INTO  FILE 


DIMENSION  LIKE (-50:50 ,0:50) , CROSS (-50:50,0:50) 
REAL  LIKE 
CHARACTER*6  ANT 
CHARACTER* 10  NAM 

OPEN  (UN IT- 10 , FILE -NAM , STATUS- ' NEW' ,RECL«10000) 

WRITE  (10,700)  ANT , BW , BWINC , MAXPT 

700  FORMAT (A6 , 2F3 .3,14) 

DO  10  I -- (MAXPT-1 ), MAXPT-1 
10  WRITE  (10,701)  (  (LIKE  (I  ,J)  ,.J  =0  ,MAXPT-1 )  ) 

DO  20  I=- (MAXPT-1) ,MAXPT-1 
20  WRITE  (10,701)  ( (CROSS  ( I  ,J )  ,J»0 , MAXPT-1 ) ) 

701  FORMAT ( ( 100F8 . 3 ) ) 

CLOSE (UNIT-10) 

RETURN 

END 
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0001 

0002 

0002 

0004 

0005 

0006 

0007 

0006 

0009 

0010 


FUNCTION  D1ST (X,Y) 


COMPLTi  DISTANCE  FROM  CENTER  to  arbitrary  matrix  cell 


INTEGER  X,Y 

□  1ST* (((X!**2+(Y)**2)«*0.5) 

RETURN 

END 


00CL  FUNCTION  GETVAK (QUERY) 

0002  : 

0003  •*>* 


0004 

c 

RETRIEVE  ANSWER  FROM  USER 

0005 

C**»* 

UUOo 

% 

w 

0007 

CHARACTER* (*)  QUERY 

0008 

2 

0009 

WRITE (6,700) 

0010 

730 

FORMAT (X) 

0011 

2 

0012 

WRITE (6 , * )  QUERY 

0013 

READ (5,*)  GETVAR 

0014 

2 

0015 

RETURN 

0016 

END 

# 


0001 

LOGICAL  FUNCTION  YSSNO (QUERY) 

0002 

c 

Gu63 

<2****1 

#  *  *  * 

0004 

THIS  FUNCTION  RETURNS  TRUE  OR  FALSE  TO  A  YES/NO  QUESTION 

* 

0005 

*  *  *  • 

0006 

c 

0007 

CHARACTER* (*)  QUERY 

0003 

CHARACTER* 1  INPUT 

0009 

2 

0010 

WRITS (6 , 700 ) 

0011 

700 

FORMAT (X) 

0012 

n 

0013 

WRITE (6,*)  QUERY 

0014 

READ (5,791)  INPUT 

0015 

70  L 

FORMAT (Al) 

0016 

c 

0017 

I F ( INPUT . EQ . ’ Y  1  )  YESNO= .TRUE . 

0013 

IF ( INPUT . EQ . 1 y ' )  YESNO= . TRUE . 

0019 

C 

0020 

RETURN 

0021 

END 
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0001 

0002 

0  JO  3 

0004 

0005 

0006 

000' 

0003 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0053 

0059 

0060 


PROGRAM  SYSEFF 


RECEIVER: 

RXDIST  -  RX  ANTENNA  GROUND  DISTANCE  TO  TARGET  GROUND 
RXAZIM  -  RX  AZIMUTH  ANGLE 

TRANSMITTER: 

TXMAXRANGE  -  GROUND  DISTANCE  FROM  FRONT  OF  MATRIX  TO  REAR 
TXMID  -  GROUND  DISTANCE  FROM  BEGINNING  Or'  MATRIX  TO  TARGET 

BOTH: 

NUMBER  OF  MATRIX  CELLS  IN  GROUND  MATRIX 
RELATIVE  SIZE  OF  GROUND  CELL  IN  DIRECTION 
ORTHOGONAL  TO  BEAM 
ANTENNA  NAME 

ANTENNA  INCLINATION  ANGLE 
BEAM  INCIDENT  ANGLE 
ANTENNA  HEIGHT 

ANTENNA  TO  TARGET  SLANT  RANGE 
ANTENNA  GROUND  DISTANCE  TO  FRONT  OF  MATRIX 
ANTENNA  BEAMWIDTH 
ANTENNA  INCREMENT  ANGLE 
NUMBER  OF  DATA  POINTS  ACROSS  BEAM 
LIKE  POLARIZATION  PATTERN  ACROSS  BEAM 
CROSS  POLARIZATION  PATTERN 
VERTICAL  POLARIZATION  FLAG 
HORIZONTAL  POLARIZATION  FLAG 


SIZE: 

YSIZE: 

TX (RX) ANT: 

TX (RX) THETAO : 

TX (RX) NOMTHETA: 
TX(RX) HEIGHT: 

TX ( RX ) NOMRANGE : 
TX (RX) LOC: 

TX (RX)BW: 

TX (RX) BWINC: 

TX (RX)MAXPT: 

TX (RX ) LIKE: 
TX(RX) CROSS: 

TX  (RX)VERPOL: 
TX(RX)HORPOL: 


1 

2 

1 

2 


COMMON  /LARGE/SIZE, YSIZE, AREA 

COMMON  /TX/TXANT ,TXTHETA0 ,TXNOMTHETA , TXHEIGHT ,TXNOMRANGE , 
TXMAXRANGE, TXMID, TXLOC  ,TXBW ,TXBWINC , TXMAXPT , 

TXLIKE , TXCR05S , TXVERPOL , TXHORPOL 
COMMON  /RX/RXANT , RXTHETAO , RXNOMTHETA , RXHEIGHT , RXNOMRAHGE , 

RXLOC , RXAZIM , RXBW, RXBWINC , RXMAXPT , 

RXLIKE , RXCROSS , RXVERPOL , RXHORPOL 
COMMON  /TERRA IN /SI GW , S IGVH  , S IGHV , SIGHH , ANGVV  ,  ANGVH , ANGHV , ANGHH 


DIMENSION  TXLIKE (-50:50,0:50)  ,TXCROSS (-50:50,0:50) 

DIMENSION  RXLIKE (-50: 50, 0:50)  , RXCROSS (-50 :  50 ,0 : 50 ) 

REAL  SCATVV (0:8) ,SCATVH(0:3)  ,SCATHV(0:3) ,SCATHH(0:3) 

REAL  PHASVV (0:3)  ,PHASVH (0:8)  ,PHASHV(0:8) ,PHASHH(0:3) 

CHARACTER* 6  TXANT , RXANT , TERNAM 

CHARACTER* 8  RUNNAM 

CHARACTER*11  RUNFILNAM 

CHARACTER*10  TXFILNAM, RXFILNAM 

CHARACTER* 1  MONOSTAT IC , RXNOMOMATCH 

INTEGER  TXVERPOL, TXHORPOL, RXVERPOL, RXHORPOL 

LOGICAL  PRINT 


RESET  PRINT  FLAG  FOR  OUTPUT  TITLES 
PRINT*, FALSE. 


SYSSFF 

0061 

READ <5,702) RUNNA” 

0062 

"32 

FORMAT (A8) 

0063 

c 

0064 

2 

GET  DATA  FROM  RON  SET  FILE 

0065 

2, 

0066 

OPEN (UNIT-20 , FILE -RUN NAM// ' . IN ' .STATUS- ’OLD 1 ) 

0067 

READ  (20,200) MONOSTATIC , RXMOMOMATCH , TXANT , RXANT  , 

006-: 

i 

TXMID.TXHSIGHT, 

006y 

2 

RXSLANT , RXNOMTHETAMIN , RXN3MTHETAMAX , 

0070 

3 

RXTHET AOMIN , RXTHETaOMAX , 

0071 

4 

RXAZIM, 

0072 

5 

SIZE, YSIZE, TERN  A>1, 

0073 

6 

TXVERPOL , TXHORPOL , RXVERPOL , RXHORPOL 

0074 

200 

FORMAT (2A1, 2A6.10F1 4. 7,A6,4I1) 

0075 

c 

0076 

CLOSE (UNIT=20) 

0077 

2 

0078 

2 

GET  TERRAIN  FILE 

0079 

2 

0080 

OPEN  (UNIT-20 ,FILE»TERNAM// ' .TER  1 , STATUS- ’OLO' ) 

0081 

z 

0082 

z 

OPEN  OUTPUT  FILE 

0033 

w 

0084 

OPEN  (UNIT-30, FILE-RUNNAM//' .OUT' .STATUS- ' NEW ) 

0085 

z 

0086 

z 

IF  SYSTEM  GEOMETRY  IS  MONOSTATIC 

0087 

z 

GET  ISOTROPIC  SCATTERING  COEFFICIENTS,  PHASE  ANGLES 

0083 

2 

0089 

IF  (MONOSTATIC. EQ. ’Y* .OR. MONOSTATIC. EQ. ’y' )  THEN 

0090 

2 

0091 

READ  (20,201)  SIGVV ,SIGVH ,SIGHV ,SIGHH , 

0092 

1 

ANGVV , ANGVH , ANGHV , ANGHH 

0093 

201 

FORMAT  (8F14.7) 

0094 

2 

0095 

2 

OUTPUT  TERRAIN  CHARACTERISTICS 

0096 

z 

0097 

WRITE  (30,300)  S I  GW,  ANGVV,  SIGVH,  ANGVH, 

0093 

1 

SIGHV , ANGHV ,SIGHH , ANGHH 

0099 

300 

FORMAT  ('1','  TERRAIN  CHARACTERISTICS’,//, 

0100 

1 

’  W  -  SCATTERING  COEFFICIENT  =  ’,F8.3, 

0101 

2 

’  PHASE-  ' ,F8 . 3 ,/, 

0102 

3 

’  VH  -  SCATTERING  COEFFICIENT  -  ’,F8.3, 

0103 

4 

•  PHASE-  ' ,F3 .3 ,/, 

0104 

5 

'  HV  -  SCATTERING  COEFFICIENT  =  ',F8.3, 

0105 

5 

’  PHASE-  • , F3 . 3 , / , 

0106 

7 

'  HH  -  SCATTERING  COEFFICIENT  =  ',F3.3, 

0107 

8 

’  PHASE-  ' ,F8 . 3 ) 

0108 

C 

0109 

2 

CONVERT  SCATTERING  CHARACTERISTICS  FROM  DB 

0110 

z 

0111 

SIGVV-10** (SIGVV/10 ) 

0112 

SIGVH-10** (5IGVH/10 ) 

0113 

SIGHV-10** (3IGHV/10 ) 

0114 

SIGHH-10** (SIGHH/10) 
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SYSEFF 

0115 

- 

0116 

2 

0117 

3 

0113 

0119 

- 

0120 

0121 

202 

0122 

2 

3123 

0124 

2 

0125 

0126 

150 

0127 

0123 

203 

0129 

c 

0130 

w 

0131 

s* 

0132 

0133 

304 

0134 

0135 

0136 

0137 

0138 

0139 

w 

0140 

r* 

0141 

c 

0142 

10 

0143 

0144 

n 

0145 

C 

0146 

0147 

0148 

0149 

c 

0150 

C 

0151 

C 

0152 

C 

0153 

0154 

w 

0155 

C 

0156 

0157 

0158 

0159 

500 

0160 

0161 

0162 

c 

0163 

“2 

0164 

0165 

0166 

0167 

GET  DATA  FOR  B I STATIC  CAS 5 
ELSE 

READ  (23,202)  AZIM.TXISC 
FORMAT  (2Fo . 3 ) 

GET  SCATTERING  COEFFICIENTS,  PHASE  ANGLES 
DO  150  I >0 , 3 

READ  (20,203)  SCATVV  ( I ) ,3CATVH ( I ) .SCATHV ( I ) ,SCATHH ( I ) , 
PHASVV(I) , PHASVH ( I ) ,PHASHV(I) ,PHASHH(I) 
FORMAT  (8F14.7) 

OUTPUT  BISTATIC  TERRAIN  CHARACTERISTICS 
WRITE  (30,304) 

FORMAT  ('1',’  TERRAIN  CHARACTERISTICS' ,//, 

'  BISTATIC  TERRAIN  -  VARIES  WITH  INCIDENT  ANGLE') 

ENDIF 

CLOSE  (UN IT»20 ) 

RETRIEVE  TX  ANTENNA  INFO 
TXFI LN AM-TXANT// ' . ANT ' 

CALL  GETANT (TXANT ,TXFILNAM ,TXBW,TX8WINC ,TXMAXPT ,TXL IKE , TXC ROSS ) 
RETRIEVE  RX  ANTENNA  INFO 
RXFILNAM-RXANT//' .ANT' 

CALL  GETANT ( RXANT , RXFI LN AM , RXBW, RXBWINC , RXMAXPT , RXLIKE , RXCROSS ) 

BEGIN  CYCLING  THROUGH  POINTS  TO  BE  CALCULATED 
START  WITH  RX  ANTENNA  INCLINATION 

DO  90  RXTHETA0»RXTHETAOMIN , RXTHETA0MAX ,10 

TELL  THE  OUTSIDE  WORLD  WHERE  YOU  ARE 

OPEN  (UNIT«50 ,FILE«RUNNAM// ' .PRG ' , STATUS” ' NEW ' ) 

WRITE  (50,500) RXAZIM,RXTHETAO 

FORMAT  ('  JUST  STARTED  AZIMUTH  «  ',F10.3,'  RXTHETA  =  ',F10.3) 
CLOSE  (UNIT*50) 

CHECK  IF  THETAO  AND  NOMINAL  THETA  ARE  TO  BE  THE  SAME 

IF  ( RXNOMOMATCH . EQ . ' Y ' . OR . RXNOMOMATCH . EQ . ' v ' >  THEN 
RXNOMTHETA*RXTHETAO 
GOTO  50 


ENDIF 
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3  1 3  L 
j  L  6  £ 
Oi.0  4 


017  3 
3  L  74 
0 17  3 
0176 
01“’ " 
0 1  7  ' 

0179 

0130 
0131 
0182 
0133 
013-5 
0135 
0136 
0137 
0183 
0139 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0  L99 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0203 
0209 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
02L7 
0219 
0219 
0  220 


cycle: 

DO  9  J 


CALC J LA  r-  R \ 


THROES:; 

SXNOMTHII 

/ERTICAL, 


beam  incident  angles 

TA=RXNOMTHETAMIN ,  KXNQMTHE  TrtMAT  ,  1  - 
HORIZONTAL  DISTANCES  FROM  TARGE  I 


53  RXNOMRANGE=RX5LAN f 

RXHE IGHT=RXNOMRANGS *COSD ( RXNOMTHETA ; 

RXDI5T=RXNOMRANG5*3  INJ  (RXNOMTHETA  i 

G  CHECK  TO  SEE  IF  TX  AND  RX  ARE  IN  SAME  LOCATION 

C  IF  So,  SET  TX  LOCATION  EQUAL  TO  RX  LOCATION 

I"  (MONOSTATIC .EQ. ' Y ' .OR. MONOSTAT I C . EQ . 1  y  ' i  THEN 
TXTHETAO=RXTHETA0 
TXNOMTHET A “RXNOMTHETA 
TXNOMRANGE=RXNOMRANGE 
TXMID “TXNOM RANGE 
TXMAXRANGE=2*TXMID 
TXHE IGHT  =  RXHE I GHT 
RXLOC=TXMID-RXD 1ST 
TXLOC=RXLOC 

ELSE 

TXMAXRANGE=2*TXMID 

TXNOMTHETA= ATAND (TXMI D/TXHE IGHT ) 

TXTHSTA0 “TXNOMTHETA 

TXNOMRANGE=TXHE IGHT /COSO (TXNOMTHETA) 

TXLOC=0 . 0 

RXLOC=TXMID— RXDIST 

C  GET  B I STATIC  TERRAIN  AT  THIS  INCIDENT  ANGLE 

S  IGVV  =  10  *  * (SCATVV ( RXNOMTHETA/IO )  /10) 
ANGVV*PHASW  ( RXNOMTHETA /10  ) 

3  IGVH  =  10  *  * (SCATVH ( RXNOM  THETA/1 0 ) / 1 0  ) 
AHGVH»PHASVH (RXNOMTHETA/IO) 

SIGHV=10*» (SCATHV (RXNOMTHETA/10) /10) 
ANGHV=PHASHV (RXNOMTHETA/10 ) 

S  IGHH  =  10  *  * (3CATHH (RXNOMTHETA/10) /10) 
ANGHH=PHASHH (RXNOMTHETA/10) 

END  IF 

30  IF  ( .NOT. PRINT)  THEN 

C  OUTPUT  TRANSMITTER  INFO 

C 

WRITE (30 , 301 ) TXANT , TXVERPOL  ,  TXHORPOL , TXNOMRANGE 
301  FORMAT!///,'  TRANSMITTER  DATA',//, 


1 

'  ANTENNA 

TYPE 

> 

,  6  X  ,  A  6  , 

2 

'  ANTENNA 

VERTICAL  POLARIZATION 

> 

,nx,i: 

3 

'  ANTENNA 

HORIZONTAL  POLARIZATION 

> 

,  L1X  ,  1 1 

4 

'  ANTENNA- 

->TARGST  RANGE 

> 

,F12. 3, 
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SYSEFF 

0221 

z 

0222 

2 

OUTPUT  RECEIVER  INFORMATION 

022i 

2 

3224 

WRITS (30,302) RXANT , RXVERPOL , RXHORPO 

L , RXNOMRANGE 

0225 

332 

FORMAT (///,'  RECEIVER  DATA’,//, 

0226 

1 

'  ANTENNA  TYPE 

> ' , 6 X , A6 , / , 

0227 

2 

'  ANTENNA  VERTICAL  POLARISATION 

>'  ,11X,U,/, 

0228 

3 

'  ANTENNA  HORIZONTAL  POLARIZATION 

>',11X,I1,/, 

0228 

4 

'  ANTSNNA->TARGET  RANGE 

> ' »F12. 3 ,/) 

0230 

z 

0231 

** 

OUTPUT  DATA  HEADER 

0232 

0233 

WRITS  (30,303) 

0234 

^33 

FORMAT  Cl1 ,2X, 'TRANSMITTER ' ,17X, 'RECEIVER' , 

27X, 'POWER  RETURN’,/, 

0235 

1 

IX, 14 (’  +  ')  ,  10X  ,  22  ('  +  ')  , 10X , 37 ('  +  ') ,/ , 

0236 

2 

'  INCLIN ' ,3X, 'INC ID' ,10X, '  AZ',3X,'IN 

CLIN' , 3X , ' INC ID ' ,10X, 

0237 

3 

'  VV ' ,8X, 'VH ' , 3X , ' HV ' , 8X , ' HH ' ) 

0238 

z 

0239 

PRINT-. TRUE. 

0240 

ENDIF 

0241 

2 

0242 

90 

CALL 

RUNONE 

0243 

z 

0244 

CLOSE 

(UNIT-30) 

0245 

CLOSE 

(UNIT-90) 

0246 

z 

0247 

END 
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SUBROUTINE  RUMON£ 


THIS  SUBROUTINE  RUNS  ONE  SNAPSHOT  0?  THE  OEPOLARI ZATION 
EFFECTS  FOR  A  PARTICULAR  RX  AZIMUTH  AND  INCIDENCE  ANCLE 


LIST  OF  VARIABLES: 


TX (RX ) GROUND : 


TX ( RX ) THETATAB : 
TX (RX) PHITA3: 

TX ( RX ) THETAPTAB : 
TX (RX ) PHI PTA3 : 
TX(RX) POLAR: 


(1)  GROUND  MATRIX  HOLDING  LIKE 
POLARIZATION  ANTENNA  PATTERN 

(2)  GROUND  MATRIX  HOLDING  CROSS 
POLARIZATION  ANTENNA  PATTERN 
THETA  ANGLE  OF  BEAM  ON  CELL 
PHI  ANGLE  OF  SEAM  ON  CELL 

THETA  PRIME  ANGLE  OF  BEAM  CELL 
PHI  PRIME  ANGLE  0?  3EAM  CELL 

(1)  COSINE  (PSI)  OF  BEAM  ON  CELL 

(2)  SINE  (PSI)  OF  BEAM  ON  CELL 


1 

2 

1 

2 


COMMON  /LARGE/SIZE,YSIZE,AREA 

COMMON  /TX/TXANT , TXTHETAO ,TXNOMTHETA .TXHEIGH ", TXNOMRANGE , 
TXMAXRANGE , TXMID  ,-TXLOC  , TXBW ,TXBWI NC , TXMAXPT , 

TXLI KE ,TXCROSS .TXVERPOL , TXHORPOL 
COMMON  /RX/RXANT , RXTHETAO , RXNOMTHETA , RXHEIGHT , RXNOMRANGE , 

RXLOC , RX A  Z I M , RXBW , RXBWI NC , RXMAXPT , 

RXLIKE , RXCROSS , RXVERPOL , RXHORPOL 
COMMON  /TERRA IN/SIGVV ,  SIGVH , SIGHV , SIGHH , ANGVV , ANGVH , ANGHV , ANGHH 


CHARACTER*6  TXANT , RXANT , RUNNAM 

DIMENSION  TXLIXE(-50:50,0:50) ,TXCROSS (-50:50,0:50) 
DIMENSION  RXLIXE (-50 : 50 ,0 : 50 ) .RXCROSS (-50:50,0:50) 

DIMENSION  TXGROUND(200, -200:200, 2) , RXGROUND  (200,-200:200,2) 
DIMENSION  TXTHETATAB(200, -200:200) .TXPHITAB  (200 ,-200 : 200 ) 
DIMENSION  RXTHETATAB (200,-200: 200) , RXPHITA3  ( 200 , -200 : 200 ) 
DIMENSION  TXTHETAPTAB(200, -200:200) .TXPHIPTAB (200 ,-200  :  200 ) 
DIMENSION  RXTHETAPTAB (200,-200:200)  , RXPHIPTAB (200 , -200  :  200 ) 
DIMENSION  TXPOLAR (200,-200:200,2) .TXRANGE ( 200  ,-200 : 200 ) 
DIMENSION  RXPOLAR (200,-200:200,2) , RX RANGE (200,-200:200) 
INTEGER  STATUS, TXMAXPT, RXMAXPT 
INTEGER  TXVERPOL.TXHORPOL, RXVERPOL, RXHORPOL 


CALL  PAINT (TXGROUND ,TXLOC .TXLIKE , TXCROSS , 

1  TXMAXPT , TXBW , TXBWI NC , TXMI D , TXMAXRANGE , 

2  TXNOMTHETA, TXTHETAO, TXHEIGHT,  TXRANGE, 

3  TXTHETATAB, TXPHITAB, TXTHETAPTAB, TXPHIPTAB, ) 


C  DEPOLARIZE  TX  COMPONENTS 


CALL  SURDEPOL ( TXTHETAO , TXTHETATAB .TXPHITAB , 

1  TXTHETAPTAB , TXPHIPTAB , TXPOLAR ) 

'OLD  ’GROUND'  OUT  OVER  ENTIRE  MATRIX 

CALL  DOUBLE (TXGROUND .TXTHETATAB .TXPHITAB .TXRANGE .TXPOLAR) 
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•PAINT'  SURFACE  WITH  RECEIVER 

CALL  PAINT ( RXG ROUND , RXLGC , RXLIKE < RXCSOSS , 

rxmaxp; ,axau ,  rxbwinc  ,txmiu  ,txmaxrange, 

RXNOMTHC r A ,  RXTHETAO  ,  RXHEIGHT  ,  RXRANGE  , 

RXTHETA  TAB , RXPH I TA3 , RXTHETAPTA3 , RXPH I ?  r  A3 ) 

DEPOLARIZE  RX  COMPONENTS 

CALL  SURDEPOL (RXTHETAO ,  RXTHETATAB , RXPHITAB , 

RXTHETAPTA3 , RXPHI PTA3 , RXPOLAR ) 

UNFOLD  RXGROUND  ALONG  LINE  OF  SYMMETRY 

CALL  DOUBLE (RXGROUND , RXTHETATA3 , RXPH ITA3 , RXRAHGE , RXPOLAR ) 

ROTATE  RX  POSITION  ACCORDING  TO  RXAZIM 

IF  (RXAZIM. NE. 0.0)  THEN 

CALL  ROTATE (RXAZIM, RXGROUND, RXRANGE, 

RXTHETATAB , RXPHITAB , RXPOLAR ) 

ENDIF 

INTEGRATE  OVER  SURFACE 

CALL  INTEGRATE  (TXGROUND ,TXPOLAR,TXTHETATAB , TXr  HITAB ,TXRANGE  , 
RXGROUND , RXPOLAR , RXTHETATAB . RXJHITAB , RXRANGE , 
TXNOMRANGE,TXVERPOL,TXHORPGl , 

RXNOMRAMGS , RXVERPOL , RXHORPOL , 

WTOTPOW ,  VHTOT  POW ,  HVTOTPOrt ,  HHTOT  PO  W , 

S I  GW  ,  S IGVH  ,  S IGHV ,  S IGHH  ,  ANGW  ,  ANGVH ,  ANGHV  ,  ANGHH ) 

CONVERT  OUTPUT  POWER  TO  DB 

IF  (WTOTPOW. NE. 0.0)  VVTOTPOW«10*ALOG10  (ABS  (WTOTPOW)  ) 

IF  ( VHTOT POW. NE .0.0)  VHTOTPOW«10»ALOG10 ( ABS (VHTOTPOW)  ) 

IF  ( HVTOT  POW . N  E . 0 . 0 )  HVTOTPOW*10*ALOG10 (ABS (HVTOTPOW) ) 

IF  (HHTOTPOW.NE .0 .0 )  HHTOTPOW»10*ALOG10 (ABS (HHTOTPOW) ) 

OUTPUT  RESULTS 

WRITE ( 30 , 3 00 )  TXTHETA0 .TXNOMTHETA , 

RXAZIM , RXTHETAO , RXNOMTHETA , 

WTOTPOW, VHTOT POW, HVTOTPOW, HHTOT POW 
FORMAT (2X,2(F5.1,3X) ,7X,3 (F5.1,3X) ,7X,4 (F7.3.3X) ) 

RETURN 

END 
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SUBROUTINE  GETANT (ANT  ,NAM,B..’,B.iISC,MAXPT  ,L 


THIS  SUBROUTINE  INPUTS  INDIVIDUAL  ANTENNA  PATTERNS  AND 
INFORMATION  ABOUT  THE  ANTENNA  FROM  THE  ANTENNA  DATA  FILE 


DIMENSION  LIKE (-50:50,0:50) , CROSS (-50:50,0:50) 
REAL  LIKE 
CHARACTERS  ANT 
CHARACTER* 10  NAM 

OPEN  (UNIT=10 , FILE “NAM , STATUS® ' OLD ' ) 

C  GET  ANTENNA  PARAMETERS 

READ (10, 700)  ANT,BW, 3WINC,MAXPT 

700  FORMAT (A6,2F3.3,I4) 

DO  10  I— (MAXPT-1)  ,MAXPT-1 
10  READ (10,701) (LIKE (1,0) , J=0 , MAXPT-1 ) 

DD  20  I— (MAXPT-1)  , MAXPT-1 
20  READ (10,701) (CROSS(I,J) , J-0 , MAXPT-1 ) 

701  FORMAT (100F8. 3) 

CLOSE (UNIT>10) 

RETURN 

END 
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SUBROUT  1  MS  PAINT  (GROUN  0  ,  LOG  ,  LI  KE ,  GROSS  , MAXPT  ,  iV.I , 

SWING,  MID,  MAXRNG, NOMTHETA, THSTAO,  HEIGHT,  RANGE, 
THETATAo , PH  I  f A3 , THETAPT Aa , PH  I PTA3 > 


THIS  SUBROUTINE  TRANSLATES  AN  ANTENNA  PATTERN  FRO.-;  ANTENNA  TO 
GROUND  AND  DETERMINES  TtfS  LOCATION  OF  EACH  GROUND  CELL  WITHIN 
THE  ANTENNA  PATTERN 


LIST  OF  VARIABLES : 


NORMXANG 

NORMYANG 


MAXYANG 

NOMTHETA 

INCANG 

XGRCELL 

YGRCELL 

XDIS 


: ANTENNA  INCLINATION  ANGLE 
:  INCIDENT  ANGLE  OF  ANTENNA  TO  GROUND  CELL 
:  ’ X ’  DISTANCE  OF  GROUND  CELL 
:'Y’  DISTANCE  of  ground  cell 

:  ’ X ’  DISTANCE  FROM  BEGINNING  OF  MATRIX  TO  CELL 


Y0AMG  : INCIDENT  ANGLE  WITH  RESPECT  TO  ANTENNA  INCLINATION 


COMMON/LARGE/S I ZE , YS I ZE , ARE A 

DIMENSION  LIKE (-50: 50, 0:50)  ,CROSS (-50 : 50 ,0 : SO ) 

DIMENSION  GROUND (SIZE , -SIZEsSI ZE , 2 ) , RANGE (SI ZE ,-S I ZE :SI ZE ) 
DIMENSION  THETATAB (SI ZE , -SI ZE:SIZE ) , PHITAB (51 ZE , -SI ZE :S I ZE ) 
DIMENSION  THETAPT Ad (SI ZE  > -SI ZE :SI ZE) , PHIPTAB (SI ZE , -SI ZE : S I ZE ) 
INTEGER  TI,TJ, MAXPT 

REAL  NORMXANG, NORMYANG, MAX YANG, NOMTHETA, LIKE 
REAL  MAXRNG ,MID , MAT , INCANG  ,  I VAL , I JVAL , LOC , LASTGA 

C  DETERMINE  GROUND  CELL  SIZE  IN  X  DIRECTION 

XGRCELL* (MAXRNG/S I ZE ) 

YGRCELL-XGRCELL/YSI ZE 
AREA>XGRCELL* YGRCELL 

W 

C  TRANSLATE  PATTERN  TO  GROUND 

n 

DO  95  I ”1 ,SI ZE 

C  FIND  INCIDENT  ANGLE 

XDIS”I "XGRCELL 
C 

C  CHECK  TO  SEE  IF  CELL  IS  BEHIND  ANTENNA 

XDIS=XDI3-L0C 

INCANG” AT AND (XD IS/HE IGHT ) 

YOANG” I NCANG-NOMTHETA 

C  DETERMINE  IF  YOANG  IS  WITHIN  BEAMWIDTH 

IF  ( ABS (YOANG *2 ) .GT . BW)  GOTO  95 
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pa::;: 

0  0  5c; 

- 

0  0  5  v’ 

DO  9 j  0  =  0, SIZE 

0  0  5  0 

00  51 

FIND  SLANT  CROSS  -.A-IGE 

0  0  'J  0 

2 

0  0  6  0 

ZD  1 5  =  J  *ZGRC£LL 

0054 

RAD= (XDIS**2rXDIS**2) ** 0 .5 

00  6  5 

SLANT* (HEIGHT**  2  i-RAD*  *  2 ' **0.5 

006o 

z 

0  06' 

** 

CALCULATE  PHI  OF  THIS  POIN: 

0  06  3 

2 

0060 

IF  (RAJ.EU.0.0)  THEN 

0070 

PHI  =0 

0071 

ELSE 

0072 

PH  I = ACOSO (XDIS/RAD) 

0073 

END  I  ? 

0074 

z 

0075 

2 

CALCULATE  THETA  OF  THIS  POINT 

0076 

2 

0077 

THETA=ACOSD (HEIGHT/SLANT ) 

0073 

2 

0079 

Z 

ROTATE  TO  ANTENNA  COORDINATE  FRAME 

0080 

2 

0081 

CALL  SURANTROT (THETA, PHI , THET A0 , THETAP , PH  I P ) 

0032 

2 

0083 

2 

IS  THETAP  LARGER  THAN  BEAMrfIDTH 

0084 

2 

0085 

IF  ( ( (90 -THETAP) *  *2  +  PH IP**2) **0.5 .GT . BW/2 )  GO 

TO  95 

0086 

2 

0037 

z 

INTERPOLATE  INTO  LIKE  AND  CROSS  TABLES  TO  GET 

GROUND  VALUES 

0Ub8 

z 

0089 

CALL  ANTINTERP (THETAP , PH I P , S I ZZ , BWINC , MAXPT , L 

IKE, CROSS, GA,CR) 

0090 

z 

009  L 

z 

FILL  GROUND  CELL  a’ITH  APPROPRIATE  VALUES 

0092 

2 

0093 

GROUND (I, J,1)*GA 

0094 

GROUND ( I , J , 2 ) =CR 

0095 

THETATAB ( I , J ) =THETA 

0096 

PHITAB (I ,  J  )  =PH I 

0097 

THSTAPTAB ( I , J ) =THETAP 

0098 

PHIPTAB ( I , J ) =PHIP 

0099 

RANGE ( I , J ) “SLANT 

0L00 

2 

0101 

90 

CONTINUE 

0102 

9  5 

CONTINUE 

0103 

2 

0104 

999 

RETURN 

0105 

END 
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ANT  INTER? 

0053  VAL1*  (LIKE  (ITH  +  INC , IPK  +  1 ) -LIKE  (IT.H,  I?H-rl)  ) 

005 L  +FTH+LIKE  ITH,IPH+i) 

006  )  GA*  (VAL1-VAL0  )  +  FPHWALO 

0061  : 

0062  T 

0063  VALO- (CROSS (ITH  + INC, IPH) -CROSS  (ITH,IPH) ) 

0064  1  +FTH+CR0S3 (ITH, IPH) 

0065  VAL1* (CROSS (ITH+ INC, I PH+1 (-CROSS (ITH,IPH+1) ) 

0065  i  *FTH+CR03S (ITH,IPH+1) 

0067  CR* (VAL1-7AL0) +FPH+VAL0 

0063  C 

0069  C 

0070  C  RETURN 

0071  END 


3  'J  B  R  3 T  I N L  3URDEP0L  (  THE  I A  J  ,  THE  I’  A  ,  PH  I  ,  THETA?  ,  ?H  I  ?  ,  POLAR 


THIS  SUBROUTINE  USES  CROC'.  -'  .tGTSNSA  ANGLES  TO  DETERMINE 
TRANSLATION  MATRIX  (COS  (Pol,,  SIS  (P31  .•  . 

(TICuLAR  GROUND  CELL 


COMMON  LARGE/S I ZS , Y3 I ZE ,  AREA 

DIMENSION  THETA (SIZE, -SIZS:3IZE,  , PH  I  ( S I ZE , -S I ZE : S I ZE ) 
DIMENSION  THETA? i 3  I  ZE , -S I Z£ : S  I  ZE j  , PH  I ? (S I ZE , -S I ZE : 3  I  IE ] 
D I MENS  I ON  POLAR (3  I ZE , -S I ZE : S I ZE , 2 } 


3INTH0=3IND (THETAJ  ! 

COSTKO =COSD (THETAO  , 

DO  93  1=1, SIZE 
DO  90  J=0,SIZE 

I  IS  THERE  A  VALUE  AT  THIS  POINT 

I?  (THETA ( I , J ) . EQ .0.0)  GOTO  95 

SINTH=3IND (THETA (I , J ) ) 

C0STH=C03D (THETA (I ,J) ) 

S  INPH=SIND (PHI  (I  ,J) ) 

COSPH=COSD (PHI  (I, J)  ) 

S INTH?=SIND (THETAP ( I , J ) ) 

C0STHP=C03D (THETAP ( I , J ) ) 

S I NPHP=S IND (PHIP (I ,J) ) 

C03PHP=C0SD (PHIP (I ,J) ) 

CALCULATE  PS I 

COSPSI=COSPH  *COSPHP+S INPH*SINPHP*SI NTHO 
S I NPS I  * (1 -COS PS I**2)**0.5 

FILL  POLARIZATION  MATRIX 

POLAR ( I , J , 1 ) =COSPS I 
POLAR (I,J,2)=SINPSI 

JO  CONTINUE 

J5  CONTINUE 


RETURN 


THETATAB  ,  .  TAB, RAN 


YLL  PRE7  i  DU  3  CALCULATIONS  WORKED  0:.’  ONLY  HALF  OF  THE  ANTENNA 
PATTER!';  T  H :  G  SUBROUTINE  ASSUMES  SYMMETRY  AND  COMPLETES  THE 
ENTIRE  PATTERN  BY  MOVING  CONTENTS  OF  A  CEL_  T  '  .  CONJUGATE 
ON  THE  OTHE.-.  SI  ■  C  OF  THE  LINE  OF  SYMMETRY 


COMMON /LARGE, '3 I ZE , Y3 I ZE , AREA 
DIMENSION  GROUND (SIZE,-SIZ2:SIZE,2) 
DIMENSION  THETATA  ,  . 3  I ZE , -3  I ZE : S I ZE ) 
DIMENSION  PHITAB(3IZE,-'SIZE:SIZE) 
DIMENSION  POLAR (SIZE, -SIZE: SIZE, 2' 
DIMENSION  RANGE (SIZE, -SIZE: SIZE) 

DO  90  1=2, SIZE 
DO  9G  J  =  i  ,  3 1 ZE 

GROUND ( I , -J , 1 ) “GROUND (I ,  J  ,  1 ! 

GROUND ( I , -J , 2 ) “GROUND (I , J , 2 ) 
THETATAB (I ,-J ) “THETATAB (I , J) 

PH I TAB (I ,-J ) “PHITA3 (I , J) 

POLAR (I ,-J,lj “POLAR (I  ,J,1) 

POLAR ( I , -J , 2 ) “POLAR ( I  , J , 2 ) 

RANGE ( I , - J ) “RANGE  ( I , J ) 

CONTINUE 

RETURN 

END 
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3E3R  j'JTINE  HOT  ATE  (  AT  IM  ,GR  j  J.V.  ,  RANGE  ,  THETA  ,  PH  1  ,  POLAR 


THIS  SUBROUTINE  ROTATES  A  RECEIVER  PATTERN  FROM  ZERO  OE 
AZIMUTr.  (WHERE  INITIAL  CALCELAT  IONS  WERE  MADE: 

::  ITS  REQUIRED  AEIM’JTH 


COMMON  /LARGE/S I ZE , YS I ZS , AREA 
DIMENSION  GROUND  (31 ZE , -S I ZE : 3  I 2E , 2 ' 
DIMENSION  THETA  (SIZE, -SIZE ’.SIZE) 
DIMENSION  PHI  {SIZE, -SIZE: SIZE) 
DIMENSION  POLAR (SIZE, -SIZE: SIZE, 2) 
DIMENSION  RANGE  (S I ZE , -3  I ZE : S I ZE ) 

CALL  SPINl (AZIM, THETA) 

CALL  SPIN!  (AZIM, PHI ) 

CALL  SPIN!  (AZIM, RANGE) 

CALL  SPIN2  (AZIM, GROUND) 

CALL  3PIN2  (AZIM, POLAR) 

RETURN 

END 
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SUBROUTINE  SPIN!  I  AZIM, MAT) 


THIS  3U3ROCTINE  ROTATES  A  SINGLE  LAYER  MAT 
TWO-DIMENSIONAL  CODE.:.  I  N  A  0  r '-.A  IE  TRANSLAT 
AND  INTERPOLATION 


COMMON, LARGS  SIZE, ¥ SIZE, AREA 
DIMENSION  TEMP  i  20 « ,-200:200) 


INTEGER  XP,XPT,Y?,XINT,YINr, YINC 
REAL  MAT (SIZE, -3IZ£:SIZE) 

REAL  X,Y,XFRAC,YFRAC,VAL0,VAL1 


DO  90  X PT  =  1 , S I ZE 


X?=XPT-SIZE/2 


DO  90  YP  =  1 , 3 1 ZE 


X=XP*COSD (AZIM) -Y?*SIND (AZIM) 
X=X+3 I ZE/2 


Y=X?*3INJ (AZIM) +¥?*COSD (AZIM) 


IE (X.LE. 1.0 .OR.X.GS.SIZE.OR.ABS (V) .GE.SIZE)  THEN 
TEMP (XPT , YP ) =0 . 0 
GO  TO  90 

END  I F 


X I NT" I NT  (X) 
XFRAC=X-X1NT 


Y I NT  =  I  NT  (Y) 
YFRAC=Y-YINT 
IF  (Y.LE.0)  THEN 
YINC=-1 

ELSE 

Y I NC  =  + 1 

ENDIF 


VAL0= (MAT (XINT,YINT) -MAT (XINT , Y 1NT+Y I NC ) ) 
VAL0=VAL0*YFRAC+MAT (XI  NT , YI NT ) 

VAL1  = (MAT (XINT+1,YI NT ) -MAT (X I NT+1 , Y I NT+Y I NC ) ) 
VAL1=VALL  * YFRAC  +MA T (X INT  + 1 , Y I  NT 
TEMP (XPT , Y  P )  =  (VAL1-VAL0) *XFRAC +VAL0 


DO  95  1  =  1, SIZE 
DO  95  J  =-S I ZE , S I ZE 


MAT (I , J ) =TEMP  (I  ,J! 


CONTINUE 


RETURN 

END 


HIS  SUBROUTINE  ACCOMPLISHES  THE  SAME 
HE ICE  OPERATING  ON  TWO-LAYER  MATRICE 


COMMON /LARGE /S I 2E , YS I ZE , AREA 
SEAL  MAT (SIZE, — S I ZE : 3 1 ZE , 2 ; 
DIMENSION  TEMP(20u, -200:200, 2) 

INTEGER  X?,X?T, YP,XINT, YINT,YINI 
REAL  X,Y, XFRAC, YFRAC , VALO ,VAL1 


DO  9/  XP?  =  L  ,  3 1 ZE 


XP=X?T-3IZE/2 

DO  90  YP  =  1  ,S I ZE 

X=XP*COSD (AZIM) -YP+SIND (AZIM) 

X=X+S I ZE/2 

Y=XP*3IND (AZIM) +YP*COSD (AZIM) 

IF (X.LE. I.O.OR.X. GE . S I ZE . OR . ABS (Y ) .GE.SIZE)  THEN 
TEMP (XPT,YP, 1 ) =0 .0 
TEMP (XPT,YP,2) =0.0 
GO  TO  90 

END  I  F 

XI NT “I NT (X) 

XFRAOX-XINT 

YINT  =  INT  (Y) 

YFRAC= Y— Y I  NT 
IF  (Y.LE.0)  THEN 
YINC  —  1 

ELSE 

YINC=+1 

END  IF 

VAL0= (MAT (XI NT, YINT  r 1 ) -MAT (X I NT , Y I NT  +  Y I MC , 1 )  ) 
VAL0=VAL0  * YFRAC  +MAT ( X  I  NT , Y I  NT , 1 ) 

VAL1= (MAT (XINT+1 , YINT, 1) -MAT  (XINT+1 , YI NT+Y I NC , 1 ) ) 
VAL1 =VAL1 * YFRAC* MAT (XINT+I , YINT , 1 ) 

TEMP (XPT , YP , 1 ) = (VAL1-VAL0) +XFRAC+VAL0 

VAL0= (MAT (XINT , YINT , 2) -MAT (X I NT , Y INT+Y I NC , 2) ) 

VALO *VAL0* YFRAC +MAT (XI NT, YINT, 2) 

VAL1 = (MAT (XINT+i, YINT , 2) -MAT ( X 1 NT  +  1 , Y I NT  +  Y I NC , 2 )  ) 
VAL1*VAL1 * YFRAC +  MAT (XINT+I, YINT, 2) 

TEMP (XPT , YP , 2) = (VAL1-VAL0 ) *XFRAC+VAL0 

CONTINUE 
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SPIN2 

0056 

0057 

00  95  I-2,SIZ£ 

0053 

DO  95  .J  =  -SI22,SIZE 

0059 

2 

0060 

MAT { I , J , 1 ) -TEMP  ( I , J , 1 ) 

0061 

HAT (I , J , 2 ) *?3MP (I , J , 2 ) 

0062 

2 

0063 

*5 

CONTINUE 

0064 

3065 

RETURN 

0066 

END 
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SUBROUTINE  INTEGRATE 

( TXGROONO , TXPOLAR , TXTHETATAB , TXPH I T AB , TXRANGE , 

2  RXGROUNO , RXPOLAR , RXTHETATAB , RXPHITAB , RX RANGE , 

3  TXMOMRANGE , TXVERPOL , TXHORPOL , 

4  KX'JOMRANGJ ,  RXVSRPOL , RXHURPOL , 

5  WTOTPOW  ,VHTOTPO«J ,  HVTOTPOW ,  H1ITOTPOW , 

6  3 1 GW , S IGVH , S IGHV  , S IGHH  ,  ANGVV , ANGVH , ANGHV , ANGHH ) 


THIS  SUBROUTINE  APPLIES  THE  RADA-  CROSS  SECTION  EQUATION  TO 
EACH  CELL  WITHING  THE  INTERSECTION  Of  THE  RX  ANO  TX  3EAM3  ANO 
ACCUMULATES  THE  RESULT  TO  FIND  THE  MEASURED  SCATTERING 
COEFFICIENT 


LIST  OF  VARIABLES: 

TXRSSP:  TX  ANTENNA  RESPONSE  MATRIX 

TXTRANS:  TX  TRANSLATION  MATRIX  (ANTENNA  TO  GROUND) 

RXTRANS :  RX  TRANSLATION  MATRIX  (GROUND  TO  ANTENNA) 

RXRESP:  RX  ANTENNA  RESPONSE  MATRIX 

SCAT:  THEORETICAL  SCATTERING  COEFFICIENT  MATRIX 


COMMON  /LARGE/SI ZE , YS I ZE, AREA 
DIMENSION  TXGROUND (SI ZE , -S I ZE : S I ZE , 2 ) 

DIMENSION  RXGROUND (SI ZE , -SI ZE : SI ZE, 2 ) 

DIMENS  ION  TXTHETATAB (SI ZE , -S I ZE : S I ZE )  ,TXPHI TAB (S I ZE , -S I ZE : S I ZE ) 
DIMENSION  RXTHETATAB (S I ZE , -S I ZE : S I ZE ) , RXPHITAB (S I ZE , -S I ZE : S I ZE ) 
DIMENSION  TXPOLAR (SI ZE , -SI ZE :SI ZE , 2 )  ,TX RANGE  (SI ZE , -SI ZE : 51 ZE ) 
DIMENSION  RXPOLAR (SIZE, -SIZE:SIZE, 2) , RXRANGE (SI ZE , -SI ZE :SI ZE ) 
DIMENSION  TXRESP (2,2)  .TXTRANS  (2 ,2)  ,  RXTRANS (2,2) ,RXRESP(2,2) 
DIMENSION  TOTPOW (2,2) 

COMPLEXES  SCAT (2,2) ,POW(2,2) 

INTEGER  TXVERPOL , TXHORPOL , RXVERPOL , RXHORPOL 
REAL  NORMAL 

C  CLEAR  ACCUMULATORS 

RXAREA-0.0 

C 

DO  5  1-1,2 
DO  5  J-1,2 

5  TOTPOW (I,J)-0.0 

DO  90  1-1 .SIZE 
DO  90  J— SIZE, SIZE 

C  CHECK  IF  RECEIVER  PATTERN  CONTAINS  THIS  CELL 

c 

IF  (RXGROUND ( I , J , 1 ) . EQ .0.0. OR. TXGROUND (I,J,1) .EQ.0.0)  THEN 
GOTO  90 

ELSE 


NORMALIZE  THE  CELL  AREA  BY  TAKING  THE  SQUARE  OF  THE  TOTAL 
LIKE-POL  PATTERN  ON  THE  GROUND  AND  MULTIPLYING  BY  CELL  AREA 


RXAR-  (RXGROUND  (X  ,J,1)  ‘TXGROUND  (I ,  X  ,1)  )  **2.0 


0053 
0  j54 
0060 
0061 
0062 
00  6  3 
0  0  6  4 
0065 
0066 
0067 
0063 
0064 
0070 
0071 
J072 
007  3 
0074 
0075 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0033 
0084 
0085 
OOdb 
0087 
0083 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0103 
0109 
0110 
0111 
0112 
0113 
0114 


C 

C 


RXAREA=RXARD\+AREA*RXAR 

ENDI? 

CHANGE  OCAS’ TER  TO  RECTANGULAR  (COMPLEX)  AND  PU  I'  INTO  MATRIX 

SCAT (1,1) -CMPLX l (3IGVV“0.5) *COSD (ANGVV/2; , 

1  (SIGW**0. 5)  *SIND  (ANGVV/2)  ) 

SCAT (1,2) =CM?LX ( <SIGVH**0.5) *COSD (ANGVH/2) , 

1  (3IGVH**0. 5) *3IND (ANGVH/2) ) 

SCAT (2,1 ) =CMPLX ( (S IGHV  **0 . 5 ) *COSD (AMGHV/2 )  , 

1  (SIGHV“0.5) ‘SIND  (ANGHV/2)  ) 

SCAT (2,2) -CMPLX ( (SIGHH**0.5) *COSD (ANGHH/2) , 

1  (3IGHH**0.5) *SIND  (ANGHri/2) ) 

RETREIVE  PSI  IN FORMA T I ON  FOR  TRANSLATION  MATRICES 

TXTRAN5 (1,1) =TXPOLAR (I ,J,1) 

TXTRAN5 (1,2) =TXPOLAR ( I , J , 2 ) 

TXTRAN3 (2,1) --TXTRANS (1,2) 

TXTRANS (2,2) -TXTRANS (1,1) 

RXTRAN3  IS  TRANSPOSED  TO  ACCOUNT  FOR  DIRECTION  CHANGE  (GROUND  TO 
ANTENNA) 

RXTRANS (1,1) -RXPOLAR ( I  ,  J  ,  1 ) 

RXTSAN5 (1,2) --RXPOLAR ( I , J , 2 ) 

RXTRANS (2,1) --RXTRANS (1,2) 

RXTRANS (2,2) -RXTRANS (1,1) 

RETRIEVE  ONE-WAV  VOLTAGE  PATTERNS  FOR  ANTENNA  RESPONSE  MATRICES 

TXRESP (1 , 1 ) -TXGROUND ( I , J , 1 ) ‘TXVERPOL 
TXRESP (1,2) -TXGROUND ( I , J , 2 ) *TXVERP0L 
TXRESP (2,1) -TXGROUND (I ,J, 2) *TXHORPOL 
TXRESP (2,2) -TXGROUND (I ,J ,1) *TXHORPOL 

RXRESP  IS  TRANSPOSED  TO  ACCOUNT  FOR  DIRECTION  CHANGE  (GROUND  TO 
ANTENNA) 

RXRESP (1,1) -RXGROUND ( I , J , 1 ) *RXVERPOL 
RXRESP (2 ,l)»RXGROUND (I ,J ,2) *RXVERPOL 
RXRESP (1,2) -RXGROUND (I ,J ,2) ‘RXHORPOL 
RXRESP (2 ,2) -RXGROUND (I , J ,1 ) *RXHORPOL 

MULTIPLY  TX  RESPONSE  MATRIX  BY  TX  TRANSLATION  MATRIX 

POW( 1,1) -TXRESP (1 ,1) ‘TXTRANS (1,1) +TXRE3P (1,2) *TXTRAN3 (2,1 ) 

POW(l ,2) -TXRESP (1,1) ‘TXTRANS  (1,2) +TXRESP (1,2) ‘TXTRANS (2,2) 

POW (2,1) -TXRESP (2,1) ‘TXTRANS ( 1 , 1 ) +TXRESP ( 2 , 2 ) ‘TXTRANS (2,1) 

POW (2,2) -TXRESP ( 2 , 1 ) ‘TXTRANS ( 1 , 2 ) +TXRESP ( 2 , 2 ! ‘TXTRANS ( 2 , 2 ) 

MULTIPLY  RESULT  BY  SCATTERING  MATRIX 


?OW(i,l)=?OW;l,l )  “SCAT  (  I  ,  l )  +P0WU.21  *  SCAT  (2,1  ■ 

POW (I, 2 > -POW (1,1' "SCAT (i,2)+POW(l,2) ‘SCAT (2. 2) 

POW (2,1)  -POW (2,1 '  ‘SCAT  (1,1)  +POWv2,2)  ‘SCAT  (2,1) 

POW  (2,2)  *?OW  (2,1)  ‘SCAT  (1,2'  +  POW  (2,2)  ‘SCAT  (2,2) 

MULTIPLY  RESULT  B Y  RX  TRANSLATION  MATRIX 

POW  (1,1) =POW (1,1) *  RXTRANS (1,1) +POW (1,2) ‘RXTRANS (2,1 ) 

POW  U,  2) -POW  (1 ,1) ‘RXTRANS ( 1 , 2 ) +POW (1 , 2 ' ‘RXTRANS (2,2) 

POW (2, l)=?OW (2,1) ‘RXTRANS (1 ,1) +POW(2, 2) ‘RXTRANS (2,1) 

POW (2, 2) -POW (2, 1) ‘RXTRANS (1 , 2 ) +POW ( 2 , 2 > ‘RXTRANS (2,2; 

MULTIPLY  RESULT  BY  RX  RESPONSE  MATRIX 

POW (1,1) =POW (1,1) ‘RXRESP (1,1) +POW (1,2)*RXRESP(2,1) 

POW (l,2)=POW(l,l) ‘RXRESP (1 ,2) +POW ( 1,2) ‘RXRESP (2,2) 

POW (2,1) -POW (2,1) ‘RXRESP (1 , 1 ) +POW (2 , 2 ) ‘RXRESP (2,1) 

POW (2, 2) -POW( 2,1) ‘RXRESP (1,2) +POW ( 2 , 2 ) ‘RXRESP (2,2) 

CALCULATE  RANGE  AND  AREA  WEIGHTING  FACTOR 

WEIGHT=AREA/( (TXRANGE ( I , J ) *RX RANGE ( I , J ) )  “2) 

SQUARE  RESULT  AND  WEIGHT 

00  75  K-1,2 
DO  75  KK-1,2 

TEMP-REAL (POW ( K , KK ) ‘CON JG ( POW ( K , KK ) ) ) ‘WEIGHT 
TOT POW (K , KK ) -TOTPOW ( K , KK ) +TEMP 

CONTINUE 

NORMALIZE  ACCUMULATOR  TO  NOMINAL  RANGE  AND  WEIGHTED  RX  AREA 

NORMAL* ( (TXNOMRANGE ‘RXNOMRANGE )  “2 ) /RXAREA 

WTOTPOW-  (TOTPOW  (1,1)  !  ‘NORMAL 
VHTOTPOW* (TOTPOW (1 ,2) ) ‘NORMAL 
HVTOTPOW* (TOTPOW (2,1) ) ‘NORMAL 
HHTOTPOW* (TOTPOW (2 ,2) ) ‘NORMAL 

RETURN 

END 


m 


TO’ 


SYSTEM  EFFECTS  ANALYSIS 


90 


SYSTEM  EFFECTS  ANALYSIS 


1.  INPUT  ANTENNA  DATA 

2.  INPUT  RUN  SET  DATA 

3.  INPUT  ISOTROPIC/MONOSTATIC  TERRAIN  FILE 

4.  INPUT  BISTATIC  TERRAIN  FILE 

5.  END 

ENTER  CHOICE  -  <2> 


RUN  SET  DATA  INPUT  ROUTINE 


INPUT  NAME  OF  RUN  FILE  <TESTRUN1> 

ENTER  MATRIX  SIZE  <50> 

ENTER  X: Y  RATIO  (I:?)  <2> 

**»  TRANSMITTER  »** 

ENTER  ANTENNA  IDENTIFIER:  <TEST99> 

FILE  FOUND 

ENTER  TRANSMITTER  POLARIZATION  <V,H)~  <V> 
***  RECEIVER  *** 

ENTER  ANTENNA  IDENTIFIER:  <TEST99> 

FILE  FOUND 

ENTER  RECEIVER  POLARIZATION  (VfH)~  <H> 


3 


RECEIVER  AND  TRANSMITTER  IN  SAME  LOCATION  (Y ,N) ?  <Y> 
ENTER  SLANT  DISTANCE  FROM  RX  TO  TARGET  <5000> 

ENTER  MINIMUM  RX  INCIDENT  ANGLE  <0> 

ENTER  MAXIMUM  RX  INCIDENT  ANGLE  <80> 

RX  ANTENNA  INCLINATION  MATCH  INCIDENT  ANGLE  (Y  ,N ) ?  <Y> 
ENTER  RX  AZIMUTH  ANGLE  <0> 

ENTER  TERRAIN  FILE  IDENTIFIER:  <MNTERR> 

FILE  FOUND 


SYSTEM  EFFECTS  ANALYSIS 


1.  INPUT  ANTENNA  DATA 

2.  INPUT  RUN  SET  DATA 

3.  INPUT  ISOTROPIC/MONOSTATIC  TERRAIN  FILE 

4.  INPUT  BISTATIC  TERRAIN  FILE 

5.  END 

ENTER  CHOICE  -  <2> 
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SYSTEM  EFFECTS  ANALYSIS 


1.  INPUT  ANTENNA  DATA 

2.  INPUT  RUN  SET  DATA 

3.  INPUT  ISOTROPIC/MONOSTATIC  TERRAIN  FILE 

4.  INPUT  3ISTATIC  TERRAIN  FILE 

5.  END 

ENTER  CHOICE  -  <3> 


MONOSTATIC  TERRAIN  DATA  FILE  INPUT  ROUTINE 

ENTER  TERRAIN  FILE  IDENTIFIER:  <MNTERR> 

FILE  FOUND 

REPLACE  ;Y,N) ?  <Y> 

ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  <0> 

ENTER  V7  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 

ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-I00> 

ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <2> 

ENTER  HK  PHASE  ANGLE  <0> 

$ 
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SYSTEM  EFFECTS  ANALYSIS 


I 

I 


1.  in pc:  antenna  data 

2.  INPUT  RUN  SET  DATA 

3.  INPUT  ISOTROPIC/MONOSTATIC  TERRAIN  FILE 

4.  INPUT  BISTATIC  TERRAIN  FILE 

5.  END 

ENTER  CHOICE  -  <4> 


i 

r.vl 

m 


i 

m 


BISTATIC  TERRAIN  DATA  FILE  INPUT  ROUTINE 


ENTER  TERRAIN  FILE  IDENTIFIER:  <BITERR> 
FILE  FOUND 

REPLACE (Y,N) ?  <Y> 

ENTER  RX  AZIMUTH  ANGLE  <45> 

ENTER  TX  INCIDENT  ANGLE  <45> 


INCIDENT  ANGLE*  0.00 
ENTER  W  SCATTERING  COEFFICIENT 
ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-10.00 
ENTER  W  SCATTERING  COEFFICIENT 
ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT 
ENTER  HH  PHASE  ANGLE  <0> 


(IN  DB)  <0> 

(IN  DB)  <~100> 
(IN  DB)  <-100> 
(IN  DB)  <2> 


(IN  DB)  <.625> 
(IN  DB)  <-100> 
(IN  DB)  <-100> 
(IN  DB i  < . 25> 


INCIDENT  ANGLE-20.00 
ENTER  W  SCATTERING  COEFFICIENT 
ENTER  W  PHASE  ANGLE  <Q> 

ENTER  VH  SCATTERING  COEFFICIENT 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT 
ENTER  HH  PHASE  ANGLE  <0> 


(IN  DB)  <1 . 25> 
(IN  DB)  <-l D0> 
(IN  DB)  <-100> 
(IN  DB)  <-1.5> 
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INCIDENT  ANGLE-33. 00 

ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  <1.375> 
ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  03.25> 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-40.00 

ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  <2.5> 
ENTER  VV  PHASE  ANGLE  <0 > 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <-5> 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-50.00 

ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  <3.12S> 
ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <-6.75> 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-60.00 

ENTER  VV  SCATTERING  COEFFICIENT  (IN  DB)  <3.75> 
ENTER  VV  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <-8.5> 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-70.00 

ENTER  W  SCATTERING  COEFFICIENT  (IN  DB)  <4.375> 
ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <-10.25> 
ENTER  HH  PHASE  ANGLE  <0> 

INCIDENT  ANGLE-80.00 

ENTER  W  SCATTERING  COEFFICIENT  (IN  DB )  <5> 

ENTER  W  PHASE  ANGLE  <0> 

ENTER  VH  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  VH  PHASE  ANGLE  <0> 

ENTER  HV  SCATTERING  COEFFICIENT  (IN  DB)  <-100> 
ENTER  HV  PHASE  ANGLE  <0> 

ENTER  HH  SCATTERING  COEFFICIENT  (IN  DB)  <-12> 
ENTER  HH  PHASE  ANGLE  <0> 
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TERRAIN  CHARACTERISTICS 


W  -  SCATTERING 

COEFFICIENT  = 

0.000 

PHASE= 

0 .000 

VH  -  SCATTERING 

COEFFICIENT  = 

-100.000 

PHASE* 

0.000 

HV  -  SCATTERING 

COEFFICIENT  = 

-100.000 

PHASE* 

0.000 

HH  -  SCATTERING 

COEFFICIENT  = 

2.000 

PHASE- 

0.000 

TRANSMITTER  DATA 


ANTENNA  TYPE  > 
ANTENNA  VERTICAL  POLARIZATION  > 
ANTENNA  HORIZONTAL  POLARIZATION  > 
ANTENNA- > TARGET  RANGE  > 


TEST99 

1 

0 

5000.000 


RECEIVER  DATA 


ANTENNA  TYPE  > 
ANTENNA  VERTICAL  POLARIZATION  > 
ANTENNA  HORIZONTAL  POLARIZATION  > 
ANTE NNA-> TARGET  RANGE  > 


TEST99 

0 

1 

5000.000 


TRANSMITTER 


RECEIVER 


POWER  RETURN 


+  +  +  +■  +  +  •♦• 

+++++++ 

+  +  ♦  ♦  +  ♦  +  ♦  ♦  4-  ♦♦♦  +  •♦•♦  +  •►♦  -f  + 

+++++++++ 

4-  +  +  ♦  +  +  +  +  +  +  ♦  4-  +  +  +  +  4-  +  +  ♦  +  ♦  +  +  +  +  +  ♦ 

INCLIN 

INCID 

AZ 

INCLIN 

INCID 

vv 

VH 

HV 

HH 

0.0 

0.0 

0.0 

0.0 

0.0 

0.000 

-13.443 

0.000 

0.000 

XO.O 

10.0 

0.0 

10.0 

10.0 

0.000 

-13.744 

0.000 

0.000 

20.0 

20.0 

0.0 

20.0 

20.0 

0.000 

-13.807 

0.000 

0.000 

30.0 

30.0 

0.0 

30.0 

30.0 

0.000 

-13.907 

0.000 

0.000 

40.0 

40.0 

0.0 

40.0 

40.0 

0.000 

-13.922 

0.000 

0.000 

50.0 

50.0 

0.0 

50.0 

50.0 

0.000 

-13.918 

0.000 

0.000 

60.0 

60.0 

0.0 

60.0 

60.0 

0.000 

-13.865 

0.000 

0.000 

00.0 

70.0 

0.0 

70.0 

70.0 

0.000 

-13.870 

0.000 

0.000 

80.0 

80.0 

0.0 

80.0 

80.0 

0.000 

-14.129 

0.000 

0.000 

$ 
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TERRAIN  CHARACTERISTICS 

BISTATIC  TERRAIN  -  VARIES  WITH  INCIDENT  ANGLE 


TRANSMITTER  DATA 

>  TEST99 

>  I 

>  0 

>  4999.995 


ANTENNA  TYPE 

ANTENNA  VERTICAL  POLARIZATION 
ANTENNA  HORIZONTAL  POLARIZATION 
ANTE N N A- > TARGET  RANGE 


RECEIVER  DATA 


ANTENNA  TYPE  > 
ANTENNA  VERTICAL  POLARIZATION  > 
ANTENNA  HORIZONTAL  POLARIZATION  > 
ANTENNA- > TARGET  RANGE  > 


TEST99 

0 

1 

5000.000 


TRANSMITTER  RECEIVER 


+++  +  •*■+  +  +  +  •<■  +  +  +  +  ++  +  +  +++++  +  +  +  +  +  +  +•  +  +  +  +  +  + 


INCLIN 

INCID 

AZ 

INCLIN 

INCID 

45.0 

45.0 

45.0 

0.0 

0.0 

45.0 

45.0 

45.0 

10.0 

10.0 

45.0 

45.0 

45.0 

20.0 

20.0 

45.0 

45.0 

45.0 

30.0 

30.0 

45.0 

45.0 

45.0 

40.0 

40.0 

45.0 

45.0 

45.0 

50.0 

50.0 

45.0 

45.0 

45.0 

60.0 

60.0 

45.0 

45.0 

45.0 

70.0 

70.0 

45.0 

45.0 

45.0 

80.0 

80.0 

$ 


POWER  RETURN 

++++++++++++++++++++++++++♦+++++++ 


vv 

VH 

HV 

HH 

0.000 

-12.797 

0.000 

0.000 

0.000 

-15.871 

0.000 

0.000 

0.000 

-17.215 

0.000 

0.000 

0.000 

-18.457 

0.000 

0.000 

0.000 

-19.887 

0.000 

0.000 

0.000 

-21.238 

0.000 

0.000 

0.000 

-22.396 

0.000 

0.000 

0.000 

-23.264 

0.000 

0.000 

0.000 

-24.090 

0.000 

0.000 

] 
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Ideal  Antenna  Pattern. 

Light  Line  -  Like  Polarization  Response. 
Dark  Line  -  Cross  Polarization  Response. 


rid  Antenna  Pattern. 

it  Line  -  Like  Polarization  Response. 

it  Line  -  Cross  Polarization  Response. 


Figure  C- 1 .  Measured  Like-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 

« _ * _ 4  _ iL  A®  HP _ - _ iAA. _ ¥ _ U _  OA® 


Measured  Cross-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  0°,  Transmitter  Incidence  =  30°. 


Measured  Cross-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  90°,  Transmitter  Incidence  =  30°. 


Measured  Cross-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  135°,  Transmitter  Incidence  =  30°. 


Measured  Cross-Polarized  Scattering  Coefficient  vs.  Receiver  Incident  Angle: 
Receiver  Azimuth  =  180°,  Transmitter  Incidence  =  30°. 
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